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Abstract: We compute how one-loop bulk effects renormalize both bulk and brane effec- 
tive interactions for geometries sourced by co dimension- two branes. We do so by explicitly 
integrating out spin-zero, -half and -one particles in 6-dimensional Einstein-Maxwell-Scalar 
theories compactified to 4 dimensions on a flux-stabilized 2D geometry. (Our methods apply 
equally well for D dimensions compactified to D — 2 dimensions, although our explicit for- 
mulae do not capture all divergences when D > 6.) The renormalization of bulk interactions 
£f) ■ are independent of the boundary conditions assumed at the brane locations, and reproduce 

standard heat-kernel calculations. Boundary conditions at any particular brane do affect 
£Sj , how bulk loops renormalize this brane's effective action, but not the renormalization of other 

distant branes. Although we explicitly compute our loops using a rugby ball geometry, be- 
cause we follow only UV effects our results apply more generally to any geometry containing 
co dimension- two sources with conical singularities. Our results have a variety of uses, in- 
cluding calculating the UV sensitivity of one-loop vacuum energy seen by observers localized 
on the brane. We show how these one-loop effects combine in a surprising way with bulk 
back-reaction to give the complete low-energy effective cosmological constant, and comment 
on the relevance of this calculation to proposed applications of codimension-two 6D models 
to solutions of the hierarchy and cosmological constant problems. 
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1. Introduction 

Does the vacuum have energy? If so, does it gravitate? Much of what we do not understand 
about quantum field theory is contained in these deceptively simple questions because calcu- 
lations robustly indicate the vacuum should have lots of zero-point energy, yet cosmological 
observations indicate that this energy gravitates very little. 

This disagreement is particularly sharp if there are only four dimensions because then 
the Lorentz invariance of the vacuum makes its energy equivalent to a cosmological constant, 
which acts as a homogeneous and isotropic obstruction to having the comparatively flat 
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universe in which we appear to live. By contrast localized energy sources need not curve 
all directions equally; for example General Relativity predicts the world-sheet geometry of a 
cosmic string to be flat regardless of the value of its tension, whose main effect is to curve 
the dimensions transverse to the string world sheet (and in particular to produce a conical 
singularity at the string position) [[[J. 

This observation suggests exploring whether extra dimensions can help understand how 
the vacuum energy gravitates, such as if we were to live on a four-dimensional analog of a 
cosmic string within a spacetime having a few relatively large dimensions. Fewer dimensions 
are better in this context since bulk fields fall off less quickly with distance, and so allow all 
branes to compete with the bulk regardless of how far apart they are in the extra dimensions. 1 
This has sparked the construction of several such brane- world systems, within one or two 
^j, [|, |7], ||, [| extra dimensions. The explicit solutions identified in this way have the 
property that their on-brane geometry is flat despite having large on-brane tensions. 

Of course this in itself does not provide a solution to the cosmological constant problem, 
which also requires an understanding of why these geometries should be robust against quan- 
tum corrections. In particular, although solutions with flat on-brane geometries exist, so too 
do solutions with curved on-brane geometries. What is required is an understanding of how 
the on-brane curvature depends on physical choices for the bulk and the branes, and whether 
these choices remain stable against renormalization as high-energy modes are integrated out. 

The most progress understanding these issues has been made for 6D models, for which 
fairly general yet explicit calculations can be made. Although no special magic is found for 
non-supersymmetric models |10||, theories with a bulk described by 6D supergravity have 



very attractive features |j, 11, [12[]. In particular, although explicit solutions with on-brane de 
Sitter geometries are known 2 , a sufficient condition for the absence of on-brane curvature 
is the absence of a coupling between the branes and a particular bulk field: the scalar dilaton 
that is related to the graviton by 6D supersymmetry j|, [l(|. This is attractive since it is 
the kind of condition that is stable against arbitrary loops involving only on-brane particles 

Such a condition would not be stable against bulk loops, however, since the brane must 
couple to the metric and the metric couples to the dilaton. So bulk loops must provide an 
important part of any naturalness story, particularly tracking how loops of heavy bulk states 



1 More than two extra dimensions might exist, but would not be relevant at low energies if they were much 
smaller than the ones of interest here. 

2 These solutions are interesting in their own right as a counter-example to the many no-go theorems for de 
Sitter solutions to higher-dimensional supergravity 
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contribute to the low-energy effective vacuum energy. A missing step in this story is an 
explicit calculation of the UV sensitivity of Casimir energy calculations on the 6D geometries 
of interest. (See, however, [[?]] for an assessment of bulk UV sensitivity for Ricci-flat 
geometries - including also the gravity sector, but in the absence of branes.) 

This paper and its companion [18| close part of this 'bulk gap', by computing explicitly 



the UV-sensitive part of the corrections to both brane and bulk interactions obtained by inte- 
grating out low-spin (spins zero, half and one) bulk fields in an extra-dimensional spacetime 
sourced by two codimension-two branes. This paper presents general results for arbitrary 
low-spin fields, while the companion specializes the results to the case where the bulk matter 
comes from a 6D supergravity. Our restriction to low-spin fields is a temporary one due to 
the technical complications of diagonalizing the full gravity-sector spectrum in the geometries 
of interest, and we intend to report on calculations using the full spectrum at a later date. 3 
The background geometry of the extra dimensions with which we work is a rugby ball 
[Q, ||, ||. This is described in §|2| below, and is basically a flux-stabilized sphere with conical 
singularities at both of its poles corresponding to the back-reaction of codimension-two branes 
located there. In a nutshell, §|| argues that the Casimir energy obtained by integrating out a 
bulk field of 6D mass m has the generic form 

where r is the rugby ball's radius and a is related to its defect angle, S, (see below for more 
precise definitions) by a = 1 — 8/2n. In the limit of large mr the dimensionless function, J 7 , 
becomes 



F(mr, a) ~ Tq{q) + 



(mr) 6 , , (mr^ 



6 



s-i(a) — s (a) + (mr) si(a) — s 2 (a) H 



ln(mr) , 
(1.2) 

where J~o(a) is m- independent and the s/%(a)'s are calculated explicitly for spins zero, half 
and one in §|j. These constants contain the dependence on a in this limit, and so also encode 
the dependence of the answer on the boundary conditions of the bulk fields near the branes 
situated at the poles. Our results reduce in special cases to results in the literature for spheres 



mm- 



The logarithm appearing in eq. (|1.2[) is slightly more complicated than what normally 



arises for loop calculations with tori |17], 24], and §y shows this ultimately can be traced to 



the existence of a number of effective interactions involving the curvature (or the background 



3 See also |n| for a partial calculation of the Casimir energy of the gravity sector and for the bosonic 
part of the spectrum on a rugby ball. 
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flux), that happen to vanish when evaluated for tori. These effective interactions arise both 
in the bulk and on the branes, and are renormalized by quantum loops of bulk fields. The 
resulting running produces the logarithmic coefficient, and because it can be traced to the 
renormalization of UV divergences this logarithmic running (and the power-law dependence 
on m that pre- multiplies it) captures the dominant sensitivity to very heavy bulk loops that 
we seek. 

Furthermore, in the special case that the renormalized interactions are in the bulk la- 
grangian, the coefficients S{ are known for arbitrary geometries using very general Gilkey-de 
Witt methods [25, 26 [, a result we summarize in Appendix |A|. We check that our results 
reduce to these general results in the appropriate limit: a — > 1. A well-known property of 
the Gilkey-de Witt coefficients is that their contributions to bulk counterterms never depend 
on boundary conditions. Physically this is because they capture the effects of very short- 
wavelength modes, and because these see only local properties of the fields they don't 'know' 
about conditions imposed at the boundaries. More precisely, the only UV divergences that 
directly involve the boundary conditions are those that renormalize the brane action at which 
the boundary conditions are applied. 

For the rugby ball this implies the renormalization of bulk interactions is identical to that 
for the sphere, and once these universal bulk counter-terms are subtracted we can separately 
identify how the brane action is renormalized by bulk loops. This exposes how these renor- 
malizations depend on boundary conditions, and how they contribute to the coefficients Sj. 
Again, because short- wavelength modes cannot know about conditions at distant boundaries, 
our results for the renormalization of brane-bulk interactions are not specific to rugby balls, 
and apply equally well to any codimension-2 brane situated within a 6D geometry, and give 
rise there to a conical singularity. Explicit results for the s^'s generated by loops of low-spin 
bulk fields are given in §|], where we also see that our results agree (when appropriate) with 
known divergence calculations for spacetimes with conical singularities 

Finally, §|| identifies a subtlety that brane back-reaction introduces when using standard 
Casimir energy calculations to identify the effective 4D cosmological constant as seen by a low- 
energy observer on the brane. The quantity V computed above is the standard fare of Casimir 
energy calculations: it is the (negative of the) loop-corrected 1PI effective lagrangian density 
evaluated at the background classical geometry (in this case a rugby-ball) . And in the absence 
of branes this quantity is also the effective cosmological constant seen by a 4D observer, at 
least for systems where the zeroth-order 4D geometry is flat. This is because the 1PI action of 
the low-energy 4D effective theory has the generic form, C e g/y/—g = —A + curvature terms, 
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which becomes —A if evaluated at the classical 4D background (which was assumed to be 
flat). As §[| shows, what changes in this argument once brane back-reaction is included is 
that it is no longer sufficient to evaluate the 1PI action at the uncorrected classical spacetime, 
even at first order in the loop corrections. This is most clear for systems where back-reaction 
cancels a classical brane tension, since in this case renormalizations of the tension should also 
be cancelled in the same way. 

In summary, what we present here as new is: an explicit calculation of the divergent 
part of the Casimir energy obtained from loops of low-spin bulk fields in a 6D geometry 
compactified on a flux-stabilized 2D rugby ball. We explicitly show how these divergences are 
renormalized into bulk and brane counter-terms, and compute how the corresponding effective 
interactions run as a result. Finally, we show how these renormalizations can be matched to 
the low-energy 4D effective theory, to see how they feed through to the cosmological constant 
as seen by a low-energy 4D observer. 

While we think these calculations can have a variety of applications to loop effects in 
extra-dimensional spacetimes, our main application is to use them to examine how super- 
symmetry ameliorates the UV sensitivity of the vacuum energy, as described in a companion 



paper [18]. Brief comments on the relevance to using codimension-2 branes to address the 
cosmological constant problem are summarized in §0. 



2. Bulk field theory and background solution 

We begin by summarizing the field content and dynamics of the bulk field theory of interest: 
D-dimensional matter (with spins zero, half and one) coupled to gravity. We then describe 
two-dimensional compactifications of this system in the presence of d = D — 2 dimensional 
brane sources. We do not work within the probe limit, and so explicitly include the back- 
reaction of these sources on the bulk geometry. In our explicit one-loop calculations we 
specialize to the case D = 6 and d = 4. 

Field content and action 

The fields of interest consist of a metric g MN , plus a collection of scalar fields gauge 
potentials A a M , and spin-half fermions ip r . We imagine the scalars and spin-half fields to 
transform under the gauge group, respectively represented on these by hermitian generators 
(t a )% and (T a y s . 
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The bosonic part of the classical lagrangian for these fields is, in the Einstein frame: 

4= = - J^2 R - \ QiiWDuP ^ V - K N F hMN ~ U(<f>) , (2.1) 

where D M denote gauge-covariant derivatives for the scalars, F^ N is the field strength for the 
gauge potentials, and the functions Gij{(f>), Habi^ 1 ) an d U{4>) are to be specified. 

Rugby-ball compactifications 

In general some scalars carry gauge charge and so having these be nonzero in the background 
would give some gauge bosons masses. Since in what follows our main interest is in back- 
ground configurations for the massless gauge fields (though we do include massive gauge 
fluctuations about these backgrounds), when solving the classical field equations we assume 
that all nonzero background scalars do not carry the charges of the nonzero background gauge 
fields. In this case, the equations of motion which follow from the lagrangian, eq. ( |2.1[) , are: 

1 



GMnp - lu ah M F a MN F bMN - u ti (<f>) = o, 



D M (n ab (4>)F bMN )=0, (2.2) 



Rmn + JGiMDutf D N 4? + K 2 U ah ($) F a MP F N 



2k 1 



D 



-U{4>) + \'H ab { ( t>)F a PQ F b ^ 



9mn — . 



where 



iP := g MN 



d M d N <p - v p MN d P <i>> + T-j,^) d M ^ k d N 4> 1 ] , (2.3) 



and T^ IN and 7j^(</>) respectively denote the Christoffel symbols constructed from the space- 
time metric, g M N, and the target-space metric, Gij{4>)- 



The simplest compactifications [31] are found using the Freund- Rubin ansatz [32] for 
which (jf is a constant and 

(g^vix) 0\ dF. (° ° 1 (2 4) 

V g mn (y)J ^ \0 / e mn (y) J 

where g^ v is a maximally-symmetric Lorentzian metric (i.e. de Sitter, anti-de Sitter or flat 
space), and g mn is the metric on the two-sphere, S2, whose volume form is e mn . The quantity 
/ appearing in the background gauge field — which could be any one of the gauge fields 
present in the theory — is a constant. All other fields vanish. 



4 Our metric is 'mostly plus' and we follow Weinberg's curvature conventions p9[ , which differ from those 
of MTW |3C|] only by an overall sign in the definition of the Riemann tensor. 
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The gauge potential, A m , that gives rise to such a field strength, F mn , is the potential of 
a magnetic monopole and so satisfies a flux-quantization condition in the presence of charged 
matter. We denote the background gauge coupling constant by 



1 



u 



(2.5) 



with 7i the component of H a \, corresponding to the nonzero background flux F mn . With this 
definition, requiring gauge transformations be single-valued for charged matter fields with 
charge qg for integer q gives the quantization condition 



2vriV 



F = 4nr 2 qf, 



(2.6) 



where N = 0, ±1, ... is an arbitrary integer and r is the sphere's radius, in our conventions 
satisfying R mn = —9mn/r 2 - Quantization requires the normalization constant, /, to satisfy 

N 



f 



2qr 2 



(without brane sources) 



(2.7) 



for all matter fields in the theory. If all charged fields have the same charge it is conventional 
to define g so that q = 1. When more than one nonzero charge is present we choose q = 1 for 



the smallest nonzero charge (say) and then for any second charged field with q ^ 1 eq. (2.7) 
requires there to exist another integer N such that N/q = N\ is also an integer, and so all 
charges are rational multiples of the smallest one. 

With the above ansatz — and using the identities F mp F n p = f 2 g mn and F mn F mn = 2/ 2 
— the field equations boil down to the following three conditions: 



R 



K 



D-2 

,2 



K 



D 



f 2 

2U + (D-3)^ 
9 



N 2 



2U gftu 



D-2 \iq 2 g 2 r 4 

2 

'2U + (D-S) 



and 



2U-^ 




9 1 


, i 



2U 



D 

N 2 



4q,2g2 r 4 



{21 



and so 



With / fixed by flux quantization, eq. (|2.7|) , these equations can be solved for (j)' 
also g(<j>) and U(4>) — as well as r and the curvature in the D — 2 directions spanned by the 
(/uz^) coordinates. 

Brane sources 

The solutions as outlined so far describe an extra-dimensional 2-sphere supported by flux, 
with metric 



ds 2 = r 2 (d9 2 + sin 2 9 dip : 



(2.9) 
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without the need for brane sources |31[. However a class of solutions with brane sources 
can be included very simply ||, just by allowing the angular coordinate to be periodic with 
period (p ~ ip + 2na with a not necessarily equal to unity Geometrically this corresponds 
to removing a wedge from the sphere along two lines of longitude and identifying points on 
opposite sides of the wedge [||, This introduces a conical singularity at both the north and 
south poles, with defect angle 5 = 2tc(1 — a), a geometry called the rugby ball. 5 

Physically, such a construction corresponds to the introduction of two identical brane 
sources, one situated at each of the singularities, with Einstein's equations relating the defect 
angle to the properties of the branes. Concretely, take the action of the brane to be 6 

S b = - f d d x L b 

Ab 
2ff 2 

where the ellipses denote other terms involving two or more derivatives, and the coefficients 
Tb, Ab, Bb, Cb and so on could depend on the extra-dimensional scalars <p l . However, the 
existence of a rugby ball solution does require the derivative with respect to <p of the total 
brane Lagrangian to vanish at the background, because the near-brane boundary condition 



with L b = T b - ^4 e mn F mn + ■ ■ • , (2.10) 



for the bulk scalars requires [15] 

where p denotes proper distance from the brane. 

For conical singularities, the near-brane boundary conditions for the metric imply [K], 
|j~5| l the defect angle at the brane's position is given by 

5 b = K 2 L b . (2.12) 

A rugby-ball solution requires identical branes at each pole, 7 for which 

l_ a = ^± = 4 G 6 L ± , (2.13) 

where L± is the lagrangian for either of the source branes. 

The presence of the brane sources complicates the flux quantization condition in two 
important ways. The first complication arises because the resulting defect angle changes the 



5 'Rugby ball' is used rather than 'football' to avoid a cultural ambiguity in what the shape of a football is. 
6 A more covariant way of writing the term linear in F mn is as the integral of the Hodge dual, * F, over the 
d-dimensional brane world-sheet JIo| . 

7 See [| for solutions with conical singularities that can differ at the two poles. 
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volume of the sphere, which appears in the flux-quantization condition when integrating over 
the bulk magnetic field, 

F = 4Trar 2 f. (2.14) 



'S 2 {a) 

The second complication arises because the branes themselves can carry a localized flux, given 

by 

•» = £■ (2 ' 15) 

For two identical branes the total flux localized in this way is <E>i := ^2 b &ib, in terms of which 
the flux-quantization condition becomes 

2vriV 



2vriV- 



'52(0) 

The condition on the normalization constant / is then 



VA+ / F = 2tt$i +4vrar 2 /- (2.16) 
V Js 2 (a) 



f = Q , (with brane sources) (2-17) 
2qr z 

where 

M := oj(N - $) (2.18) 



and (for later convenience) uj := 1/a, $ := (and also $6 := 



16 J- 



Control of approximations 

Since our entire discussion takes place within the semi-classical approximation we must de- 
mand all fields vary slowly enough to trust the low-energy effective-field-theory approximation 



[p3| , 34 1 for whatever (possibly a string theory) ultimately provides its ultraviolet completion. 

In practice, without knowing the details of this UV completion, we ask fields to vary 
slowly relative to the length scale, £, set by the gravitational coupling: k 2 = £ D ~ 2 . Since 
(barring unnatural cancellations) eqs. ( |2.8D imply 1/r 2 ~ k 2 V ~ k 2 N 2 /g 2 r A , we also see that 
r 2 ~ K 2 N 2 /g 2 ~l/n 2 V and sor»l also implies g 2 <C N 2 l D ~ A and V < £~ D . 

Finally, once brane sources are included we must also demand them not to curve exces- 
sively the background geometry, and for branes with tension T this requires k 2 T -C 1. For 
rugby-ball geometries this ensures the defect angle satisfies 5 <C 2ir. 



3. General features of bulk loops 

We now turn to the size of one-loop quantum fluctuations about the rugby-ball background 
just discussed. In particular, we calculate the UV-divergent part of loops computed for the 
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various fields (p l , ip r and A a M expanded about this background. Notice for these purposes 
that we need not restrict ourselves to the loop contributions of fields that are nonzero in the 
background. 

The main assumption we use to compute a field's one-loop contribution to divergences is 
to suppose that its kinetic operator can be written in the form 8 A = — □ + X + m 2 , for some 
choice of local quantity A (perhaps a curvature or background flux) and a squared mass, 
m 2 . This is sufficiently general to include most of the spin-zero, -half and -one particles of 
interest in later sections. As we shall argue in more detail below, because we restrict to the 
UV sensitive part of the calculation our results apply more generally than just to rugby-ball 
geometries. On the other hand, our calculation is insensitive to effective interactions that 
cannot be distinguished using only a rugby-ball geometry, such as differences between R 2 
and R MN R MN interactions (which are indistinguishable for spheres), or those interactions like 
V M R\7 M R involving gradients of the curvature (which vanish for spheres). 

The other main restriction to our calculations that emerges in subsequent sections is the 
need to avoid fields whose fluctuations mix nontrivially with those of the metric. These must 
be avoided because for them it is not straightforward to show that A takes the desired form. 
In particular this precludes our computing the effects of those fields that are nonzero in the 
classical background. 

Finally, for simplicity our final expressions also specialize to D = 6 and d = 4. Although 
our methods work equally well for any D and d = D — 2, our explicit evaluations only capture 
all one-loop divergences for D < 6. For D > 6 they give only a subset (the most divergent) 
of UV divergences. 

One-loop calculations 

Writing the generator of 1PI correlators as T = S + S, then the UV-sensitive part of one 
loop quantum corrections, E, can be calculated using heat-kernel methods. For a real D- 
dimensional field of mass m we have the 1-loop quantum action 9 

= d^Vi-ioop = "("If ^Tr Log ( " D " + f + , (3.1) 

where (— 1) F = +1 for bosons and —1 for fermions, A is a local combination of background 
fields (such as curvatures and background fluxes) and Vi_i 00 p denotes the effective 1-loop 
scalar potential (or Casimir energy density) in d = D — 2 dimensions. We assume that for 

8 For one-loop purposes this can also be done for fermions by working with the square of Ip. 

9 We return below to how the kinetic operator for higher-spin fields can be put into the form — □ + X + m 2 . 
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the compactification of interest the higher-dimensional d'Alembertian splits into the sum of 
d- and two-dimensional pieces: d D = + Cta- 

Anticipating the need to dimensionally regularize ultraviolet divergences we write the 
spacetime dimension as d = d — 2e, where d is an integer and e — > at the end of the 
calculation, so that 

where m 2 n denote the eigenvalues of — D2 + X in the compactified space, and we Wick rotate 
to Euclidean signature using d d k = id d k E . Using the identity In X = — J °°(ds/s) exp(— sX) 
(which is valid up to an infinite constant that is independent of X) we then have 

Vi-ioop = "("If \f £ / ^ jH y exp [s (kl + m 2 + mj n )] , (3.3) 

where [i 2 in the exponential is absorbed into a rescaling of s. Performing the d d k E integral 
using 

!^e-^ = (±X /2 (3 4) 

gives 

,2e 



Vi-ioop = -(-i) f 2(47r)d/2 e y ^1+^72 exp [ _s (™ 2 + 

2e poo jj. 

= " ; ,w/o / "TtL e-*( mr ) S(t) , (3.5) 



where t = s/r 2 is dimensionless, 

5(f):=(-l)^exp[-tA 3(1 ] , (3.6) 



and the dimensionless quantities A, n are defined by 



m 2 „ := ^ . (3.7) 



In the examples of interest 1 jr is the generic Kaluza-Klein scale for the compactification. 
In the appendices, we show that the function S(t) has the following small-t limit: 

S (t) = ^T + '^f + s ° + S V2^ + *i * + s 3/2* 3/2 + s 2 t 2 + 0(t 5 / 2 ) , (3.8) 
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where the coefficients Sj depend on the spectra, A, n , and so also on the spin of the particle in- 
volved, and in principle also on the boundary conditions used near any branes situated within 
the background geometry. Much of what follows is devoted to computing these coefficients 
explicitly for the fields and boundary conditions of interest. 
Using this small-i expansion in eq. fl3.5|) gives 



1— loop 



2- 



2 ( 4vrr 2)d/2 J Q 



dt e 



-t(mr) 2 



S-l j s ~l/2 j Sp Sj/2 
j-2+d/2 ^3/2+d/2 -/-l+d/2 + £l/2+d/2 



Si 



+ 



53/2 



+ 



S2 



fd/2 1 £-l/2+d/2 1 £-l+<f/2 



+ 0(t -2+d/ 2) 



m d /j 2£ 
"2(47r) d / 2 



(mr)- 2i r(i-d/2), 



(3.9) 



which shows that the result diverges when d ^ d for all terms, s,, for which i — d/2 is a 
non-positive integer. 

At this point we specialize to d = 4 so that all divergences are captured by Si with i < 2. 
For this case we have 



Vl-loop 



m d /i 2£ 



(mr) h - — — r 

6 2 (mr) 1 (mr) 4 



r(-d/2) + (finite as d -> 4) . (3.10) 



(47r) d / 2 

Using x 4 ~ d T(— d/2) = (4 — d)" 1 + In a; + finite, the divergent part of Vi_i oop emerges as 



Vl-loop 



(47rr 2 ) 5 



4 — d \m 



+ V 



(3.11) 



where Vt is finite and /i-independent in the limit d — > 4, and 



s_i 



C := Y'( mr ) 6 ~~ ~7j~( mr ) 4 + si(mr)' 



«2 



(3.12) 



What is important in what follows is that the coefficient C depends on the independent 
external variables that control the properties of the background geometry, like r, N, and 
a. Subsequent sections use this dependence to extract more information about the effective 
interactions that renormalize these divergences. 



Renormalization 

Ultraviolet divergences are renormalized, as usual, into counter-terms in the action, and in 
the setup of interest here this action has both bulk and brane contributions. The goal of the 
next subsections is to separate each of these types from one another. 
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Bulk counterterms 

The crucial feature of bulk counterterms that allows them to be separated from brane coun- 
terterms is their insensitivity to the boundary conditions satisfied by the bulk fields in the 
vicinity of the branes. This is most easily seen if they are computed using Gilkey-de Witt heat- 
kernel techniques [25, [2(| — see Appendix ^] — since this calculation is explicitly boundary- 
condition independent (for bulk counterterms). Physically this arises because divergences are 
sensitive only to modes with wavelengths much shorter than the physical size of the space- 
time, since this both ensures their effects are captured by local interactions and that they are 
too short to build correlations between points in the bulk and distant boundaries. 

Because they do not depend on the boundary information, bulk counterterms can be 
computed using only the bulk geometry, without making reference to brane properties. In 
particular, they may be obtained by specializing the coefficients Sj to the special case where 
the background geometry has no defect angle: a = 1. 

The interactions to be renormalized are found by writing the most general local bulk 
lagrangian consistent with the given field content and symmetries, organized into a derivative 
expansion: C B = C B Q+£ B 2+C B i • • •• Restricting to terms that are nonzero in the background, 
this expansion gives 



£-b2 

£>bA 

£b6 



-gU 
~9 



F? i ' TP rpMN 



-g 



2k 2 4g 2 
k Car 



Itr MN r H it 



C«3 R 3 + • • • 



Ch 2 "T5 2 

K 



(3.13) 
(3.14) 
(3.15) 



and so on. Here we define 



a R R + 2b R R MN R MN + c R R 



MNPQ 



j^MNPQ 



(3.16) 



with a R + b R + c R = 1 so that R = R 2 when specialized to a spherical geometry (for which 
RmnpqR mnpq = 2-R mn i? mn = R 2 = 4/r 4 ). A similar, but more elaborate, definition is used 

3 

for R . Calculations on a sphere can only track the overall renormalization £ fl 2, Cr 3 an d not 
how the separate parameters such as a R , b R and c R renormalize, but — as summarized in 
Appendix |A| — for the bulk contributions these separate renormalizations are known explicitly 
for general geometries from earlier work [^] . A similar story also holds for terms that involve 
gradients of the background scalar fields, which vanish for the rugby-ball configurations. Since 
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this includes in particular kinetic terms like (dcj)) 2 , it represents an obstruction to computing 
wave- function renormalizations for the scalars. This limits the generality of our later formulae 
for the renormalizations of scalar couplings. 

Evaluating the bulk action at the background rugby-ball solution and integrating over 
the compact directions gives 



= (4vrar 2 ) (c/-^L + — — 



(3.17) 



1 <A 



\ a c "T 



showing that [/, 1/k 2 , Ca 2 ! an d Ca 3 terms can be read off respectively from the r 2 , r°, r~ 2 , 
and r~ 4 terms in V s , while the l/g 2 and £ AH terms are identified as the coefficients of Af 2 /r 2 
and A/" 2 /r 4 , respectively. In particular, the power of M acts as a proxy for the power of /, 
and so does not appear at all if the particle in the loop does not carry the charge gauged by 
the background gauge field. 

The divergences in Vi-i 00 p are absorbed by splitting all couplings - i.e. U, 1/k 2 etc. 
- into a renormalized and counter-term part, with the infinities of Vi_i 00 p that arise in the 
bulk canceling the divergences in the bulk counter-terms as d — > 4. Once this is done the 
/^-dependence of the renormalized quantities also must cancel the explicit /x-dependence in 
the corresponding finite parts of Vi-i 00 p- And the bulk part of the divergences in Sj can be 
identified because they do not depend on the brane boundary conditions, and so are the same 
as they would have been for a calculation on a sphere. Explicitly, if we denote by s^ ph what 
the coefficients Sj would have been if evaluated on a sphere (i.e. with a = 1 and &b = 0), 
then the full result for Sj can be written 

Si -asf^ + Ssf 1 , (3.18) 

which can be regarded as the definition of 5s\ ot . (The factor of a pre- multiplying s^ ph arises 



from the integration over the extra dimensions, as in eq. (3.17).) This split is useful because 
only the divergences in the first term, as^, can be absorbed into renormalizations of the 
bulk interactions of V B , while those of 5s\ ot must be absorbed into brane interactions. 
With this definition the running of the renormalized bulk couplings satisfy 

6 . _ a / 1 \ ™4 



dU m b mh n 9/1 



_ sph,0 ( I — ____ s P h >° (1-\Qk\ 

M d/x~ 6(4vr)3 S - 1 ' ^/xUV" 2(4vr)3 S ° ' {6A ^ } 



dfj.\ k J 4(4vr) 3 1 ^ d[i 8(4vr) ; 
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and so on. Here the quantities 10 s^' , Sq 3 * 1 ' , s^ ph '° and s| p '° are computed below by 
explicitly summing over the KK spectrum on a sphere, giving results that agree with those 
obtained in Appendix [A| using general heat-kernel methods. We note at this point that since 
we derive the running of the couplings by assuming they cancel only the explicit ln(///m) 
dependence of Vi_i 00 p, our explicit formulae exclude the case where additional ^-dependence 
enters through the appearance of renormalized couplings and fields pre-multiplying the pole 
in lj{d — 4) (though the formulae are easily generalized to include this more general case). 

The renormalization of the gauge- field terms, 1/g 2 and C, AR , is similarly done by keeping 
track of those divergences involving /. This can be done, for example by comparing loops of 
particles that couple directly to the background flux with those that do not. The result is 

.. 8 (1 \ _ 2m2 ,s P h,2_ 8q 2 m 2 sph , 2 

M 0A*U 2 / (47r)3r*/ a 1 " (4tt)W 2 1 ' 1 ' 

and 

9 { k( ar \_ 2 „s P h,2_ 8g 2 aph[2 

^d^\ ~g 2 )~ (47r)3r 4 / 2 2 ~ (4vr)W 2 2 [ ' 

where the particle in the loop has charge qg, and s^ ph ' 2 represents that part of s^ ph that is 
proportional to J\f 2 . 

Brane counterterms 

A similar reasoning can be applied to brane-localized interactions, which (unlike the bulk 
counterterms) can depend on the boundary conditions used near the brane but should be 
independent of those boundary conditions imposed on distant branes. (For earlier treatments 
of divergences in the presence of conical singularities, see j28|.) 

To identify the brane contributions we first subtract the contributions of the universal 
bulk counterterms found above, using eq. ( 3.1S| ), and use 5s\ ot = ^&^ s «(6) to extract how 



each individual brane-localized interaction renormalizes. This can be done as before, by 
distinguishing those contributions that come from couplings to the background gauge field 
from those that do not. 

An additional complication in the case of brane counterterms is the necessity to disentan- 
gle which contributions come from which branes. This complication arises because although 
the KK spectra encountered below depend explicitly the flux, <£b, localized on each brane, 
they only depend on the brane tensions through the common defect angle a = 1/u) = 1— 5/2tv. 



10 Notationally, for sf ' the "sph" emphasizes that these quantities are evaluated on the sphere, and the 
superscript l k' denotes terms involving k powers of the gauge-field normalization, N '. 
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We write 



5s] 



tot 



+ 



Z^ 1 

b 



Mb) 



(3.23) 



where <5s^ ame receives equal contributions from each brane and where the Ssf^ depends on 
the explicit $ b that in general differ on each brane. Since we are interested in tracking the 
renormalization of each brane separately, the quantity of interest is 



i{b) 



+ 6s 



diff 



(3.24) 



2 ' 'Mb) 

since this is the one for which we can write 5s\ ot = 5s^ . To avoid any extra notational 
clutter, we choose to drop the subscript (b) when writing in what follows. 

Writing the most general local brane lagrangian in a derivative expansion: C b = C b Q + 
C b i + + Cb3 + • • • , and dropping terms that vanish when evaluated at the rugby ball 
background, we have 

£6o = -V=7T 6 (3.25) 

•Ab ran t-i 
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Re mn F„ 
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i7,2 



tt r M M-f 



(3.26) 
(3.27) 
(3.28) 
(3.29) 



and so on, where 7^ := g MN d ll x M d v x N is the induced metric on the brane and the e mn F mn 
terms arise covariantly as the integral of the Hodge dual, *F, over the brane world-volume. 
Evaluating these at the background solution gives a contribution from each brane of size 
, Abf 2( R b <Abf 2 ?J<ARbf 4C fl2fe n 2 ( ARb f 2 

1 ^ ~o ~o o A 



k 2g 



f 
AbN 



%Rb 



2 g 2 r 2 



+ 



nQ Ab N 2 nC, ARb N 4C fl2fc 

T~ Q ~Q A ~T~ ~ O A 1 A 



H 2 C A RbM 2 



1qg 2 r 2 Kr' 2 ' &q 2 g 2 r A 



qg 2 r A 



j 2 g 2 r & 



+ 



(3.30) 



to the 1PI potential, with the sum over branes giving V tot := V b . Requiring the counter- 
term parts of the sum of the two brane actions to cancel the remaining divergences implies 
the renormalized quantities satisfy 
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(3.31) 



- 16 - 



where both branes are assumed to be identical, Ssi is the quantity defined by eqs. fl3,18| ). 



( |3,23| ) and ( 3.24 ), for which the superscripts '0', '1' and '2' respectively correspond (as above) 
to the terms independent of, linear in, or quadratic in the background gauge field (or its 
proxy, A/"). 

What remains is to compute the coefficients s_i through S2 for the fields of interest. This 
is the aim of the next sections. 

KK mode sums 

This section now sketches how the coefficients Sj are computed, by performing the sum defining 
S(t) using a hypothetical eigenvalue spectrum, m jn , that is general enough to include most 
of the special cases of practical interest. 

We can obtain our later results by performing the sum over the mode labels, n and 
j, in either of two different ways. The most reliable (and more cumbersome) method first 
performs a Poisson resummation, which has the advantage of casting the sums in a way that 
converges more quickly for small t. The second (and easier) method avoids the complications 
of Poisson resummation, instead using zeta-function regularization to regularize the part of 
the mode sum that is non-singular as t — )■ 0. We present both methods because although 
the zeta-function technique is much simpler to use, its validity ultimately relies on the more 
complicated calculation based on Poisson resummation. More details on the equivalence of 
these two techniques can be found in Appendix pi 

To explain these two techniques, consider the following expression for the KK spectrum, 
Xj n , that is general enough to include many of the cases met in later sections, 



2 



\ jn = [j + -\n + b + \ + -\n-b_\+a) -r. (3.32) 

Here n is an arbitrary integer and j = 0, 1, 2, . . . is a whole number. The quantities a, b± 
and r are real parameters that differ for different spin fields in the loop and for different 
background flux quantum number, N. For example, in the case of a complex scalar field 
that is charged under the background flux we show in the Appendix ^| that the appropriate 
choices are a = ^, b + = \N\, 6_ = and r = j, where N is the background flux quantum. 11 
Our interest is in tracking how the sum, S(t), defined using this spectrum depends on 
the geometrical quantities N and r, as well as the rugby-ball defect angle that is encoded in 
the quantity 

u = — where a = 1 , (3.33) 

a 2tt 



11 These expressions work in a particular patch for the background monopole gauge potential, with b± 
interchanged in the opposite patch. 
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and S is the rugby-ball defect angle described in §|2[ In particular, the limiting case oj = 1 
corresponds to the sphere, for which a variety of results are known |2l| , [2^ , 23, |35| for the 
Casimir energy. For instance, in this limit and with no gauge flux (N = 0) the scalar spectrum 
becomes Xj n = 1(1 + 1), where t = j + \n\. In this case the sums can be re-ordered to give 
the usual form for scalars on a sphere: 

oo oo oo oo oo £ 

EE=EE=EE- (3-34) 

n=-oo j=0 n=—oo£=\ n \ £=Q n =—£ 

Poisson resummation technique 

We first sketch the calculation that best controls the convergence of the sums at intermediate 
steps. For clarity of explanation we do not work with the most general form for the spectrum, 
but specialize to the following special case 



oo oo 



S(oj, t) = e Tt ^2 ex P [-*(? + W M + a ) 2 ] ■ ( 3 - 35 ) 

j=0 n=— oo 

The difficulty with this sum is that it converges poorly in the regime of interest: where t is 
small. 

To remedy this we use the Poisson resummation formula, which relates the sum over a 
function f(x) to the sum over its Fourier transform J~(q). That is, if 

/oo 
dx f{x) e*** , (3.36) 
-oo 

then 

oo oo 
n=-oo k=— oo 

To apply this in the case of interest define 

fj{x) := exp [-t(j + u\x\ + a) 2 ] (3.38) 
and Poisson resum the n-sum, giving S(co, t) as 

oo oo 

S(u,t) = e Tt Y^ J2 Fj&TTk) (3.39) 

j=0 k=—oo 



where _ 

^ / ^ i pa 



a 2 e mi 1 - erf (J + zg) ) + c.c. 



(3.40) 



and 

q 



j:=Vt(j + a), ?:=t-^. (3-41) 

2 ojyt 
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Now comes the important observation, that is justified in some detail in Appendix [B|. 
Because of the factor e - ^ and the inverse power of ojy/t appearing in q, all of the terms with 
k ^ in the sum are 'regular', in the sense that their sum vanishes in the limit that either 
uj or t vanishes. Only the k = term contributes to the singular part of the small- i limit of 
S(u, t). That is, 

S(oj, t) = S sing (u, t) + S rcs (u, t) , (3.42) 



with 



OO T l I OO i- / x 

S si ^(co,t) := e ^^(0) = f-JjE 1 " erf ®) + cx 



3=0 ' 3=0 

and 



(3.43) 



S reg (u;,t) := 2e T *^^J- j (2 7 rA;). (3.44) 
fc=l j=o 

This last equality uses Fj{—q) = Fj(q), which follows from fj(—x) = fj(x). 

Appendix [B| evaluates the remaining sums explicitly in the small-i limit, giving results 
that agree with the somewhat simpler techniques we now describe. 

A simpler zeta-function method 

A somewhat simpler way to compute the small-i limit of S(u>,t) is to start with eqs. ( |3.42D 
and eq. ( |3,43| ), but not to perform the Poisson resummation for S reg (u,t). The regular part 
is instead computed by zeta-function regularizing the initial sum over n and j. 

To see how this works use the Euler-Maclaurin formula, which states that for any analytic 
function f(x), 

°° roc d 

£/»(?) =/ dxf n (x)-Y,-^ft l) (0), (3.45) 
j=o Jo i=i l - 

where Bi denote the Bernoulli numbers (of which the relevant ones are B\ = — ^, B2 = |, 
£>3 = 0, -B4 = — and B5 = 0) and where / ( - i_1 ^(0) denotes the (i — l)-th derivative of f(x) 
with respect to x, evaluated at x = 0. 

Applying this formula to the j-sums in S(ui,t), requires using the function 

fn(x) = exp [— t(x + a n ) 2 ] with a n = urn + a , (3.46) 

and so 

OO ^ ; OO _ 

S n (co,t) := = 2 V 7 i 1 ~ erf K^)) - E -f/ (4_1) (0) , (3-47) 

j=0 V i=l ' 
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where the first few fn ^(0) terms are 



/(0) (0 ) = e -K 
/< 3 >(0) = (12 t 2 a n - 8 t 3 a 3 n )e 
and /^(O) ~ C(i 3 ) for all i > 5. 



-2ta n e 



-tar, 



-tat 



/W(0) = (12 - 48 iX + 16tX)e 



/( 2 )(0) = (-2t + 4i 2 a 2 )e 



-f«r, 



3„2 



Using this in 



5 reg (w,t) 



Tt 



<S (a;,i) + 2^«S n (u;,t) 

n=l 



-fa,, 



(3.48) 



(3.49) 



and Taylor expanding the error function gives a series expression for S Icg (uj,t) that involves 
divergent sums of the form Y^n=i nk with k a non-negative integer. Remarkably, defining 
these as k), where Cr( s ) is Riemann's zeta-function, gives a finite expression, which 
agrees with S reg (u;,t) as computed using Poisson resummation. 



To this must be added S smg (uj,t), computed using eq. (|3.43|) . Once this is done the 
resulting small-i expansion for S(cj,t) can be compared with previous calculations of the 
small-t limit using Gilkey-de Witt heat-kernel expansions, when these are known. They are 
known in particular for the sphere, where co = 1, and Appendix [A| shows that they agree in 
this limit. 



4. Results for low-spin bulk fields 

We next collect results for the coefficients Si of the small-i limit for bulk fields with spins 
zero, half and one. 



4.1 Scalars 

Consider first the simplest case of a single minimally coupled real scalar field, satisfying 
0(f) = 0, that is coupled to the background gauge field with monopole number N and brane- 
localized fluxes In this case the scalar spectrum (as derived in Appendix ^| in the case 
where the north patch 12 of the gauge potential is used) is given by 



jn 



7 H re - H — n - iV + <E>_ + - 

2 1 +i 2 1 i 2 



1+A/" 2 



(4.1) 



where the superscript 's' is meant to emphasize that this (and later formulae in this section) 
applies only for scalars. The Casimir sum becomes 

oo oo f / i \ 2" 

cxp 



j=0 n=— oo 



-t[ 3 + -|n-$+| + -\n- iV + $_| + - 



• (4-2) 



2 As Appendix |c| also shows, an equivalent result is obtained if the south patch is instead used. 
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Using the results of Appendix and its notation 

F b :=\$ b \(l-\$ b \), F^:=J2K, fW:=F, G(x) := (1 - x)(l - 2x) , (4.3) 



we find the following small-t coefficients: 
f 



s s (lo,N, $ 6 ) 
s s 2 (oj,N,$ b ) 



LO 
1 

LO 

1 

LO 

1 

LO 



1 J\f 2 LO 2 

+ —(1 - 3F) 

180 24 18 v ' 



12 



- W ^$ 6G (|$ 6 |) + ^(1-15F( 2 )) 

6 



1 

504 

A i 



1 1 A' 2 / 1 A' 2 \ , ,,, 
"W + ( 90 "US.' (1 - 3F)W 



Co' 



3 A" 



24 



£>&G(|<I> 6 



+ ^-^(1-15^)) 
360 V ' 



LO 



5 A" 



120 



J2®bG(\$ b \)(l + 3F b ) 



+ i 



1 



i?(2) i?(3) 



When w = 1 and <!>>, 



1260 120 60 

0, these become 



„ s ph 



1 



„ s P h 
S 



sph, 



1 

3' 
4 
315 



„sph, 



1 

15 



sph, 2 



and s 



sph, 2 



_ _M 2 
' l ' " 24 ' 
_A^ 
40 



(4.4) 
(4.5) 

(4.6) 



(4.7) 



(4.8) 



in agreement with the result in []23| , as well as with the Gilkey-de Witt coefficients as computed 
for a 6D scalar on a sphere in Appendix ffl, using the general results found in |26| . If the 
scalar couples to the background field with strength qg, its contribution to the running of the 
leading bulk counterterms therefore is 



dU 
dfx 



m 



<> 



a 



d ( C R 2 



6(4tt) 3 

,2 



1 



m 



d_ (\_ 



60(4tt) 3 
2 q 2 m 2 
~~ 3(4^) 3 



dfi \k 2 



m 



6(4vr) 3 ' 



1 



d f k(ar 



630(4tt) 3 ' 
2q 2 



(4.9) 



5(4vr) 



It is straightforward to check that the above expression for the loop component of the renor- 
malization of the gauge coupling agrees with the result obtained directly by evaluating the 
Feynman graphs for the vacuum polarization in 6D flat space, following standard methods 
[p7|1) despite its initially unfamiliar sign. 
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Returning to general uj and to identify the brane-localized renormalizations we must 
first subtract the contribution of the bulk counterterms to obtain 5s{. As described earlier, 
because these counterterms do not depend on the boundary conditions at the branes for 
uj 7^ 1 their contribution to V contains only the trivial proportionality to 1/u = a due to the 
volume, 47rar 2 , that appears when V is computed by integrating the counterterms over the 
rugby ball. This is consistent with the overall factor of 1/u) that pre- multiplies all of the Sj 
in eqs. pL4 - 4.7 1 . We then identify 5si as prescribed by the discussion surrounding eq. ( p. 24 ) 
and find that Ss-i = 0, 
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cj — 1, and the last, approximate, equalities give the leading dependence in the 
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Here 5uj 

limit where 5uj -C 1 and |$&| <C 1. The corresponding contributions to the running of the 
brane counterterms are 
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The appearance in these expressions of non-analytic terms involving |<&&| bears some 
comment. As we see in Figure |], the signs of these terms are such that they represent 
maxima of the potential at <3?b = 0, which are cuspy in that derivatives with respect to § b 
are discontinuous at this point. We believe this discontinuous derivative arises because of 
a level crossing of the ground state at these points, for the following reasons. On the one 
hand, the explicit calculations given above reveal Vi_i 00 p as a polynomial in |3>&|, which must 
therefore grow without bound, for large |^&|« On the other hand, Vi— loop 

should be periodic 

under the replacements <&6 + 1 an d N — > N + 1, as can be seen from the invariance 

of the KK spectrum, eq. ( [4.1[ ), under this shift. (When shifting <I> + — > <I> + + 1 both N and 
n must be shifted by unity, but the shift in n is lost in Vi_i 00 p once the KK mode sum is 
performed.) This shows that the energy is not minimized for the same value of N as $5 is 
varied to sufficiently large values; instead a new vacuum with TV" — > N + 1 is energetically 
preferred once <3?& become larger than unity. Cuspy maxima in the potential can arise at the 
points where this crossover between vacua occurs. 



4.2 Spin-half fermions 



As shown in Appendix |C], the KK spectrum for a fermion that is charged under the U (1) 
whose flux supports the background rugby ball (using the north patch of the gauge potential) 
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Figure 1: Schematic plot of the potential as a function of $5 
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where ni/2 = n — a/2 and where a S {±1} specifies the 4D helicity of each component of the 
spinor, of which there are 2 (4) each in the case of a 6D Weyl (Dirac) spinor. By identifying 
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(and so A/fa- := uj(Nf a — &f a ) = A" — <r as well), we can relate these fermionic spectra to the 
scalar spectrum considered previously: 
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As discussed in Appendix |C], the corresponding small-i coefficients take different forms 
depending on whether or not |$fe| < $q. When \$b\ 5; ^0 (f° r both ^fe' s ); the mode sum over 
the above spectrum yields the following smail-i coefficients for a 6D Weyl spinor: s_-, = — 4/w, 
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(Note that the above expressions would be the ones valid when considering the limit <J>^ — 5- 
while holding oj fixed at some value 7^ 1.) When |<l?&| > <3?q, we find that s i _ 1 is unchanged, 
but that 
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where 



p b := sgn($ 6 ) = $ 6 /|$ 6 | , p := and $ 6 := w($ fe - p 6 *£) . 
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(We also use tilded versions of the notational contractions, such as Ff, := |$b|(l — \&b\), in 
the same way as is done in the previous section on scalars.) As a check of these expressions, 
we can evaluate them when |<l>f,| = $q (or simply $ = in the second case), and find they 
each give the same result. 
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These also agree with the results found using Gilkey-de Witt methods in Appendix [A|. The 
special case where q also vanishes then gives exactly one-half the result found in |23| for a 
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fermion on a sphere, as is appropriate due to our use here of 6D Weyl (rather than Dirac) 
fermions. 

For a fermion with charge qg, the corresponding contributions to the running of the bulk 
couplings are 
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Subtracting the contribution of these bulk counterterms leaves the contributions to the Gilkey- 
de Witt coefficients that renormalize the brane action. When |<&&| < <3?q, we find that <5s_i = 
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where, as before, 5u> = u — 1 and the approximate equalities give the leading result for 
5u),&b *C 1. The brane counterterms therefore renormalize as follows: 
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When instead |3>b| > we fi n d that ( as always) 5s -\ = 5s\ = and 
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Therefore, when |3>&| > $q, the brane counterterms renormalize as follows: 



/i 



^6 

3// 



d_ /A 

^ dp\g 2 



p—(^ 

djJL \ K 



A 

m 


(5oj 


2(4tt) 2 Co 
qm 2 


V 3 
Pb(l - 


(4vr) 2 w 


6 


qm 2 1 




- (4tt) 2 ^ 


3 


m 2 


/ 5a; 



+ 2F b 



m 4 / 5a; 



(4tt) 2 V 6 



+ 1*6 



+ 



2(4vr) 2 o; I 15 90 
m 2 ( Sou | <E»6, | 



5oo 2 2 5a; 3 5a; 4 F 2 

+ -4- 



u 2 F h 



45 



/i 



2(4tt) 2 
Q 



dp\g 2 J (4vr) 2 ^ 



15 3 I ' 

p 6 AA(l-a; 4 ) ATl> fe 



30 



90 



9 



(4.53) 



(4.54) 



(4.55) 



+ — 
10 



UJ 



w<i fc 



G(|* 6 |) 



a; 



30 



2q 



p b 5a; + <l b ) , 



dp " 4(4vr) 2 a; I 42 



5a; 



15(4tt) 2 

19 Slo 2 5a; 3 23 5a; 4 5w 5 



(4.56) 



1260 



F 2 F 3 uj 2 F 2 
30 15 6 



105 1200 

30 



_ 6oj6 
105 ~ 630 





f 5a; 


\*b\) 


) ~ 4(4vr) 2 


, 42 ~ 


30 ) 



(4.57) 



- 28 - 



4.3 Gauge fields 

We next state the results for the Casimir coefficient for a gauge field, provided this gauge 
field is not the field whose flux stabilizes the background 2D geometry. We consider in turn 
the cases where the 6D gauge field is massless or massive (in the 6D sense). 



Massless 6D gauge fields 

We begin with the massless case. Picking an appropriate gauge (such as light-cone gauge) 
allows the 6D gauge field to be decomposed into four components, 13 each with a spectrum 
(when evaluated using the north patch of the gauge potential) given by 
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where £ G {0, 0, +1, —1} for each of the four components. (For more details see Appendix [C] 
and |39[| .) From this, we see that two modes have the exact same spectrum as scalars (i.e. the 
£ = modes), and two modes have almost the same spectrum as scalars (i.e. the £ = ±1 
modes). Similar to before, we can identify 
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(and so Ng£ '■= oo(N g ^ — $gf^) = M + 2£) to relate these modes to the scalar spectrum: 
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However, since the resulting spectrum for the £ = ±1 modes ends up being very similar to 
that of the scalars, we shall just write 
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3 Ghosts also do not contribute to the one- loop result in this gauge. 
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so long as \<& b \ < <J>Q f . As it turns out, the | ^-dependent terms seen here in Asf { serve 
to exactly cancel any corresponding terms in 4s| in eq. ( |4.62j ) that are odd in |$&| (for the 
£ = ±1 modes only). 

Even though only S2 contributes to the Casimir energy for massless fields, we nonetheless 
also follow s—i, sq and s\ since these are useful as intermediate steps when assembling the 
contributions of a massive vector field. The corresponding bulk quantities therefore are 
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and after these are subtracted the brane renormalizations are obtained from 
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along with 5s _i = <5sf = (as usual). 

Because the renormalizations coming from Sk are proportional to m 4 ~ 2k , where m is the 
6D mass, for massless fields we need only follow the contributions of S2, ensuring the only 
nonzero renormalizations are 
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on the brane. 



Massive 6D gauge fields 

Let us now turn to massive gauge-field fluctuations, corresponding to those gauge directions 
that acquire mass because of the nonzero (but constant) values taken by some of the scalar 
fields. By assumption, these gauge fields are vanishing in the background, both because this 
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would require a more complicated ansatz than assumed here for the rugby-ball backgrounds 
pCfl , and because it would complicate the diagonalization of the metric and gauge-field fluc- 
tuations. 



In this case the linearized theory simplifies [37] if we choose light-cone gauge, as described 
in more detail in Appendix |C]. The result is that a massive gauge field leads to the 4D spectrum 
of a massless gauge field, given in eq. ( 4.59 ), plus that of a scalar provided earlier. It follows 
that the Sj coefficient of a massive gauge field are 
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where the sf are the corresponding quantities for a 6D scalar, those given in eqs. fl4.4| )-( p~7"| ), 
and we use ./V G {0, ±1} to ensure stability (see Appendix H). 

Let us now give the contribution of the massive gauge field to the running of the bulk 
couplings. Since in the sphere limit we have 
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The running of the brane couplings is similarly obtained by computing the 5si coefficients: 
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5. The 4D vacuum energy 
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The previous sections show how to compute the 1PI potential, Vi_i oop = Voo+V/, obtained by 
integrating out low-spin bulk fields, and how these divergences are renormalized into various 
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bulk and brane interactions, V B and Vb, so that V := V B + J2b + V/ is finite. In this section 
we compute the implication of the above renormalizations for the effective 4D cosmological 
constant, A, and on-brane curvature as seen by a low-energy 4D observer. In the 4D theory A 
would be obtained as the value of the low-energy 4D effective potential, V, after minimizing 
over any light scalar fields in the 4D effective theory. 

If no branes had been present, a standard result for the low-energy potential would have 
been V = V B + V/, so it may come as a surprise that once branes are included the potential 
V is not simply given by V := V, suitably renormalized. Instead we must also recompute 
the classical contribution to the low-energy cosmological constant coming from integrating 
out KK bulk modes, keeping track of how the bulk back-reacts to the renormalization of the 



source branes, along the lines of refs. [10|. Neglecting this back-reaction would be inconsistent, 



since for codimension-two systems it is known to be of the same size as the direct effects of 



the changes to the brane lagrangians themselves [|5|, 11, 12]. Indeed, it is this back-reaction 
that allows flat solutions to exist at all at the classical level, despite the large classical positive 
tensions carried by each brane. 

The logic to determining the cosmological constant of the low-energy effective theory is 
to compute within the 6D theory how perturbations to brane and bulk interactions change 
the predicted value for the curvature, R^ u , along the brane directions, and then to ask what 
cosmological constant in the effective 4D theory would give this same curvature. This is a 
special case of a 'matching' calculation between the effective theory and its UV completion 
(|35| , It can be carried out fairly explicitly for small changes about a known background 
solution, as we do below following refs. flO|| . 

An alternative route to the same end is compute the value taken by the loop-corrected 
(1PI) action, r := S + E, evaluated at the background configuration that solves the loop- 
corrected field equations, ST/Sip = 0. This lends itself well to the present purposes because the 
quantity Vi-i 00 p as computed in earlier sections is precisely such a contribution to the loop- 
corrected action evaluated at the background solution. For maximally symmetric geometries, 
the value of A in the 4D low-energy effective theory can be found by comparing the result 
for r obtained using the full 6D theory with the result for T computed using the effective 4D 
theory. 

5.1 Classical bulk back-reaction 

It is useful first to review how back-reaction works at the classical level. For the 6D theory 
we solve the bulk field equations, eqs. ( |2.2D , using the renormalized bulk couplings. The 
renormalized brane couplings enter through the near-brane boundary conditions they imply 
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for the bulk fields |15J. Our interest is in starting with a classical solution for which = 0, 
and then ask how changes as the various renormalized couplings change by small amounts. 
Of particular interest is how responds to changes of the couplings in the brane action. 

Recall that the rugby ball is sourced by identical brane actions, which to first order in a 
derivative expansion couple to the background fields by Sb = — f d 4 x^— 7 Lb with 



L b = T b --^e mn F mn + -- - , (5.1) 



2~g 2 

and where (as before) 7^ := g MN d^x M d u x N is the induced metric on the brane, and the 



ellipses denote terms involving two or more derivatives. As shown in detail in [1C], the 
nominally subdominant .A^-term can play an important role in understanding the low-energy 
effective 4D curvature in those situations [Q, || where the stabilization of the extra dimensions 
arises as a competition between brane and bulk flux. In this case the influence of the A term 
gets enhanced by the volume of the extra dimensions through its effect on the flux-quantization 
condition. 

At the classical level, the source branes back-react onto the background solutions in two 
distinct ways. First, they change the boundary conditions of the bulk fields, schematically 
relating the near-brane limit, lim p ^o(pdvp/dp), to the (appropriately renormalized) derivative 
of the brane action, 5Sb/Si/j, for any bulk field ip. In the special case that the functions T& and 
Ah defining the brane action are independent of any bulk scalar fields, then these boundary 
conditions boil down to the familiar statement that the brane induces 14 a conical defect angle 
at each brane position, of size (see eq. ( [2.12D ): 5 = n 2 Lb = k 2 {T^ — Abf /g 2 )- The second way 
back-reaction influences the background is through the flux quantization condition, eq. fl2.16|) , 
which depends on T b and A b because: (i) the defect angle changes the volume of integration 
for the flux, and (ii) because the A b term directly contributes as flux localized on the branes. 



Once both effects are included |1Q| , the flux-quantization condition generalizes from eq. (2.7) 
to 2irN/q = (Aira r 2 f + J2b A) , as before. 

What must be done is to track how the bulk solutions react to loop-induced changes, 
<5Tfc, 5Ab (and others), to see how the changes in the bulk solutions appear in the low-energy 
effective theory. 

A straw man 

Before doing so, it is worth first putting to rest a common misconception that can confuse 
issues at this point. 



14 More generally, if Tt or Ab depend on bulk scalars, back-reaction leads to a bulk curvature singularity at 
the brane positions Il5l llCj] . 
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Schematically, our goal is to compute the loop-corrected action, r = S + eS, evaluated 
at the loop-corrected background field configuration, tp = ipo + e5ip, with e being the small 
loop-counting parameter. We do so in order to compare the result when performed in 6D and 
in 4D. Working to order e gives 

^ S[ip + e5i>] + eS[^ ] + 0{e 2 ) 

' ss 



~ S[ipo] + e 



# + S[^o] +C(e 2 ), (5.2) 



and it is tempting to argue that the first of the 0(e) terms vanishes because tpQ satisfies the 
classical field equations: (5S/5ip)^ = 0. If so, then F[ip] = S[tpo] + eEfV'o] + 0( f2 )- 

Although this argument is often true, if it were true in the present instance then there 
would be no need to know how bulk fields back-react — i.e. to compute Sift — in response 
to loop corrections to the action — £ — since it would suffice to evaluate both S and S at 
the uncorrected classical solution, tpQ. If so, then the renormalized potential, V, computed in 
previous sections could be directly interpreted as the effective vacuum energy. 

To see why this argument fails it is useful to examine a concrete example. To see why 
Scsfyo + e#V>] / S e fi[tpo] at 0(e), consider the case 

S e s = - [ d 4 x^g~ | Ao + 



2^4,0 



and r 



y dW ^( A+ _|), (5 .3) 



where A = A + 5A and l/n\ = 1/k\$ + 5(1/ k%), with both corrections of order e. For 
simplicity, suppose further that Ao = 0. 

With these choices the unperturbed classical background satisfies Rq = while the back- 
ground solving the full loop-corrected equations is R\ = —Ak\ q 5K. Consequently the classical 
action evaluated at the classical background is S^ff [<7o] = 0, which disagrees with its evaluation 
at the loop-corrected solution: 



S c s[go + 5g] = - J d 4 x^g = +2 J d A Xyf^ 5 A . 



(5.4) 



Formally, the reason S e s[gQ + 6g] can differ from S c s [go] is because S e s is actually pro- 
portional to the volume of spacetime for any constant nonzero curvature, and so diverges. 
Consequently precise statements must be regularized, such as by cutting the geometry off at 
a large radius. But then the action also contains a boundary, Gibbons-Hawking, term at this 
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radius, and it is this boundary term that need not be stationary when evaluated at a solu- 
tion to Einstein's equations. A similar phenomenon was noticed for the on-shell gravitational 



action in another context in ref. [41|. 
Loop-corrected back-reaction 

We now turn to how the bulk solutions respond to the presence of the loop-generated changes 
to brane and bulk interactions. Since these changes are small within the semiclassical ap- 
proximation we can analyze this by following ref. and solving the bulk field equations 
linearized about the classical background, using the brane-bulk boundary conditions to relate 
integration constants to properties of the brane action. 

There are two qualitatively different kinds of loop-generated effects. First, renormaliza- 
tions of the bulk action captured by V B can change the features of the bulk classical solution, 
such as (but not restricted to) the generation of new terms in the bulk cosmological constant. 
Second, the renormalizations of the brane action included in Vb change the bulk solutions 
in the two ways mentioned above: changing the boundary conditions and modifying the 
flux-quantization condition. 

The calculation in can be used to include the influence of the brane renormalizations, 
once generalized to the loop corrections to the brane action found earlier in this paper. We 
now briefly describe the result of such a calculation. We first examine corrections to the 
lowest-derivative terms in the action, and return to the effects of the higher-derivative terms 
in the next section. 

It suffices to work in a simple case to make the point that back-reaction is required to 
properly infer the implications of Vi_i op for the 4D curvature. We therefore assume that 
all gauge fields that are zero in the classical background remain zero at loop level. The 
calculation is performed, without loss of generality, using Gaussian normal coordinates in 
the scalar-field target space centered at the background solution, so that the kinetic terms 



are canonical at this point. Finally, to avoid the issues of ref. [42] we imagine the branes 
only to couple to bulk scalars with vanishing bulk masses, in order not to have the brane 
and bulk compete in their implications for (/>*, and thereby preclude the existence of constant 
background configurations. 15 

We then imagine the bulk and brane actions to be chosen to admit a rugby-ball solution, 
and consider the corrections to this solution due to generic corrections to the coefficients 

k 2^ k 2 + (5k 2 ; U(<f>)^U(<l>) + 5U(<l>), 5 2 (0) ^fi^+Sg 2 ^), (5.5) 



3 Both examples of [jlOl] fall into this class 
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in the bulk and 

T b -> T b + 5T b {(j)) and A b ^A b + 5A b ((f>) , (5.6) 

on the branes. We ignore corrections to the target space metric, SQij((j)), since at linear order 
these do not contribute to the effective cosmological constant. We discuss higher-derivative 
corrections to the brane action below, but the bottom line is that they are suppressed by 
powers of the gravitational scale k. 



Following the steps laid out in ref. [10], we first compute the on-brane curvature generated 
by these perturbations, by solving the linearized 6D field equations (assuming the initial rugby 
ball to be flat in the brane directions). The the effective 4D cosmological constant, A, is read 
off as that constant that would produce the same curvature using the 4D equations of motion: 
i.e. R A = -4k|A. This gives 16 



A 



irar 



5L 



sTR^ 
fSAtot 



7TOT 



Sf 

f 



tot 



9 2 



+ 



45k z 35U 5g 



6U 

~w 

7p- 



2 5? 



+ 



U 



+ ~Zr 



(5.7) 



where / is defined (as before) by F mn = f e mn , while A to t ■= J2b^b and L tot := Y^b L b 
with L b is as given in (|5.l| ). Generically the right-hand-side of eq. (|5.7j) is evaluated at the 
loop-corrected stationary point for the bulk scalar field: + Sep 1 . When this is fixed 

from the classical bulk field equations then it suffices to evaluate eq. ( |5.7| ) a t the classical 
solution, (j>Q, since keeping b(\) % would be subdominant in the loop expansion. 

A more complicated situation arises if the bulk scalars parameterize flat directions along 
which there is no scalar potential at all in the absence of the perturbed bulk-brane couplings. 
This can occur if the brane couplings break an approximate symmetry of the bulk equations 
of motion, and can be used to generate a Goldberger-Wise [45| type stabilization of the bulk 
geometry [jD!}] . In this case it is the perturbed brane-bulk couplings that stabilize these flat 
directions, with the stabilized point, <p l , satisfying 



k 2 5L 



tot J 



4ira 



+ 



U d 6f , SUj (5g 2 



Uf 



+ 



5U 



9 2 



0. 



(5.8) 



The first term here corresponds to the case studied in ref. [10] . Because we work perturbatively 
in 8(f) 1 , some care is required if 4> l labels a flat direction of the zeroth-order equations. In 
this case (as is standard for degenerate perturbation theory) b<^ % need only be small if the 
arbitrary unperturbed point is chosen near the minimum of the loop-corrected potential. 



16 The first two terms in the last equality are ref. jn]]'s result, since the corrections to 8k, SU and 8g were 
not considered. 
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Although reasonably complicated in the general case, it is clear that the above result only 
agrees with V as computed earlier in the special case that the first term, 5L to t, dominates 
eq. ([Tt])- This is not true, in particular, in the super symmetric examples that are the subject 
of a companion paper [Il8| (for which the above formulae simplify considerably). 



5.2 Higher derivative corrections on the brane 

The above discussion is restricted to a tension and a localized flux, but in general loop 
corrections to the bulk and brane actions also include higher derivative terms, with the 
dimensions made up by factors of k. Given our experience with the flux, which contributes 
at the same level as the tension despite being down one derivative, we pause to check that 
the other terms in the derivative expansion on the brane are indeed suppressed with respect 
to the tension and flux terms. The two ways in which higher derivative terms can contribute 
is by modifying the brane flux, and by modifying the matching condition that sets the defect 
angle at the branes. 

The matching condition that relates the defect angle to the on-brane Lagrangian is more 
complicated if higher derivative corrections are taken into account. The precise statement of 



the matching condition in the general case is [15] 



[e B d p B] jT = <T + -£= {^A , (5.9) 

where the metric near the brane is assumed to take the form ds 2 = dp 2 + e 2B ^d9 2 . For 
most fields, maximal symmetry in the brane directions guarantees that the right hand side of 
this matching condition only gets contributions from the variation of y 7- 7 in the action. An 
important exception to this are curvature terms on the brane, which we now explore. 
A brane with the structure 

S b = - J d*x^ (h b + ^R^j (5.10) 

has a metric variation 

i ss b Uf R ^ v _R r \_U r = _fUR + iA ri (5n) 



V— iSguv K \ 2/2 \ 4k 2 

where the last equality uses the maximal symmetry of the on-brane directions. Effectively 
this means that in eq. ([fT?]) we should replace 5L tot — > 5L to t + (Cflft-^4) / '(2k), and re-solve for 
i?4. The result is 

" k 2 H5A tot _ §L \ _ 45k*_ 35U Sg 2 ' 
irar 2 \ g 2 tot J r 2 K 2 2U 2g 2 



(5.12) 
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From this we see that the correction that is associated with a brane curvature term is of 
higher order in the expansion in n/r 2 . 

The other contribution higher derivative terms can have is to modify the brane localized 
flux. In particular, the brane localized Maxwell term has this effect. Following the regular- 
ization of the brane flux in [^] we find a divergent contribution to the gauge potential at a 
small distance, 5, from the brane. Such divergences are normal for codimension-two matching 
[ |l5| , 45 1, and the limit 5 — > is to be taken after renormalization of the brane-bulk couplings. 



The divergent near-brane form for the bulk gauge field is in this case 
A e (d) - 



2 U + kCi&M 2 

The first term vanishes as 5 — > and so can be ignored. The other correction terms are 
also small because consistency of the semiclassical approximation requires both ng 2 <C 1 
and k5~ 2 <C 1. Again, we find corrections to the cosmological constant due to the higher 
derivative corrections to the brane are subleading in k. 



6. Conclusions 

To repeat the summary of the introduction, in this paper we present an explicit calculation 
of the divergent part of the Casimir energy obtained from loops of low-spin bulk fields in a 
6D geometry compactified on a flux-stabilized 2D rugby ball. We explicitly show how the 
UV divergences can be disentangled in order to identify separately how both bulk and brane 
interactions renormalize, and we compute the renormalization-group beta functions for each 
of the corresponding couplings. 

Although our results are general, the first applications we envisage are to the cosmological 



constant problem within the framework of supersymmetric large extra dimensions [18|. There 



are several features about this paper's results that are noteworthy for such applications. 

• The first observation is that extra dimensions are not a free lunch in themselves. 
For large m, in general all positive powers of m up to m 6 appear in the bulk loops, 
Vi-ioop of rn & r 2 . This is as expected on dimensional grounds for a contribution from 
a 6D cosmological constant. Although the m 6 term cancels whenever there are equal 
numbers of bosons and fermions, in general the subdominant Vi_i 00 p oc term is 



still present. In the absence of something special (like unbroken supersymmetry [16]), 
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the UV sensitivity we find is precisely what would have been expected generically on 
dimensional grounds. 



Second, although loops of heavy particles are dangerous in the bulk, in a world where 
all Standard Model particles live on a brane nothing really requires there to be any 
heavy bulk particles. If all bulk masses were generically of order 1/r then a single bulk 
loop need not be so dangerous, since V is of generic order V ~ C/(4-7rr 2 ) 2 . In such a 
case it would be two-loop and higher contributions that would instead be worrisome, 
since these could introduce the larger brane-localized masses of the Standard Model 
back into the result. What is required is a mechanism that suppresses higher loops, and 
here again supersymmetry is likely to be relevant. 

Should a mechanism ensure that the 4D vacuum energy is indeed of order l/(47rr 2 ) 2 (as 
we believe to be the case for a bulk supergravity fll8| ), then it is the Kaluza-Klein scale, 
1/r, that would set the size of the cosmological constant. This could be acceptably 
small if r is chosen not to be far below its current upper limit, r < 45 microns, arising 



from tests of Newton's laws at short distances [43]. (Modifications of Newton's law 
from supersymmetric large extra dimensions have been studied in Notice that the 

numerical factor of 167r 2 in the vacuum energy is an important part of why these scales 
can be compatible. 

• Finally, even should the numerical size of the Casimir energy be right, its sign is also 
important. It is in this context that the discussion of §|| is most important, since it 
shows that the sign of the low-energy cosmological constant need not simply be the sign 
of V as computed in §^. 

In any event, there is no substitute for a real calculation for which we regard the results 
of this paper as a crucial first step. 
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A. Heat kernels and bulk renormalization 

In this appendix we collect for convenience the explicit expressions for the heat-kernel coef- 
ficients for manifolds without singularities and boundaries, and specialize the results to the 
case where the two extra dimensions are a 2-sphere. 

The small-t expansion for the heat-kernel representation of the one-loop vacuum energy 
described in the main text can be evaluated for a broad class of geometries in great gen- 
erality pq ]. Such generality is possible because the small-i limit physically corresponds to 
the coincidence limit of the corresponding propagator, and this does not 'know' about the 
boundary conditions and topology of the space if the coincident points are taken far from any 
boundaries or branes. 

Gilkey-de Witt coefficients 

This section collects the results for the ultraviolet-divergent parts of the one-loop action 
obtained by integrating out various kinds of particles in 6 dimensions. To this end, consider 
a collection of ./V fields, assembled into a column vector, and coupled to a background 
spacetime metric, g M N, scalars, tp 1 , and gauge fields, A a M . We suppress the gauge and/or 
Lorentz indices to which these fields couple, leading to a background-covariant derivative, 
D M , of the form 

D M ^ = d u ^ + u M f - iA a M t a ^> , (A.l) 

where uj m is the matrix-valued spin connection, and the gauge group is represented by the 
hermitian matrices t a . For real fields the t a are imaginary antisymmetric matrices, which (for 
canonically-normalized gauge bosons) include a factor of the corresponding gauge coupling, 
g a . The commutator of two such derivatives defines the matrix- valued curvature, Y MN ^ = 
[D M ,D N ]^!, which has the following form: 

Y MN = 7&mn ~ iF mn ■ (A-2) 

Here 1Z MN is the curvature built from the spin connection oj m , which is also related to the 
Riemann curvature of the background spacetime in a way which is made explicit in what 
follows. 
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Integrating out the fields at one loop often leads to the following contribution to the 
quantum action 

%E = -(-f^Ti log(-D + X + m 2 ) , (A.3) 

where (— ) F = + for bosons and — for fermions. As before, □ = g MN D M D N and the quantity 
X is some local quantity built from the background fields. The mass matrix, m 2 , are constants. 

Our interest is in that part of S for which the functional trace is ultraviolet divergent. 
To identify the divergent part we work within dimensional regularization and so continue 
the spacetime dimension to complex values, n, which are slightly displaced from the actual 
integer spacetime dimension, 6, which is of interest: n = 6 — 2e. We then follow the poles 
in E as e —7- 0, in the usual fashion. Notice that this continuation to d = 6 differs from the 
regularization used for the Casimir energy calculation, for which the limit d — > 4 was taken. 

For 6D spaces without boundaries and singularities the resulting divergent terms are 
simply characterized. They can be written as [^6| 

S ~ = ^3 H F E r ( fc - 3 + e ) / d&x ^9 tr [m e - 2k a k ] (A.4) 
2(4vr)3 ^ J 

where T(z) denotes Euler's gamma function. The divergence as e — > is contained within the 
gamma function, which has poles at non-positive integers of the form T(— r + e) = (— ) r /(r!e) + 
• • •, for e an infinitesimal and r a non-negative integer. The coefficients, a k , are known matrix- 
valued local quantities constructed from the background fields, to which we return below. The 
trace is over the matrix indices of the at , of which there are N' = N d with N counting the 
number of fields and d being the dimension of the relevant Lorentz representation. 

The above expression shows that for massless fields {m = 0) in compact spaces without 
boundaries and singularities in 6 dimensions the divergent contribution is proportional to 
tr [03], so the problem reduces to the construction of this coefficient. By contrast, for massive 
fields there are also divergences proportional to tr[m 6 ao], tr[m 4 ai], tr [m 2 a2] and tr[<Z3]. 
(Notice that the freedom to keep m 2 within or separate from X implies that the divergence 
obtained from computing just (23 using X m = X + m? gives the same result as computing ao 
through 03 using only X.) 

An algorithm for constructing the coefficients a k is known for general X and D M , involving 
and the result for the first few has been computed explicitly [gfj] and can be given as a closed 
form in terms of X, background curvatures and the generalized curvature Y MN . The first few 
coefficients are given explicitly by [26]: 17 

a Q = I (A.5) 



7 In comparing with this reference recall that our metric is 'mostly plus' and we adopt Weinberg's curvature 
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ax = --(R + 6X) (A.6) 
o 

a 2 = t^k (2RabmnR abmn — 2R MN R MN + 5R 2 — 12 OR) 

+-RX + -X 2 - -DX + — Y MN Y MN (A.7) 

6 2 6 12 17 

a 3 = 1 ^—18 \3 2 R + YJD M RD M R - 2D L R MN D L R MN - AD L R MN D N R ML 
+9D K R MNLP D K R MNLP + 28RDR - 8R MN OR MN + 24R M N D L D N R ML 



I ion |— 1 riMNLP p3 i p p pA/JV P P TjABMN 

+ 1ZH MNLP \-JK — It -\ — — rCrC MN K —Itlt ABMN rC 



,,„ _ 35 „, 11 ,. ,. „„„ 14 

*-mnlp>— >R 9 ~ 3 2 

l pAf p piVL ^4 pA/JV pKI p I ^ pAf p pi 1 

"1 — n jv n-ML -ft -ft fimknl i ^~ -ft at i^mklp Jl 

_ ^4 pAB p rjMNKL ^0 r>A M p r>BKNP 

"g - -ft A/iV f^ABKL -ft "g" -ft S JV ftAif AfP -ft 

+ 3^o ( 8£> *' y ™ flMyM + 2£)M ^m + i2y MJV ny MJV (A. 

1 OVM V-^ V K R TtMNKL \r V J- A P? M V VJVif 

— 1^2 jv 2 K 2 jf-Oil 1 MiV f XL i~ *-ft N I MK 1 

-5RY MN Y MN - 6D 2 X + + 3(XD M Jf D M X - 60X 3 

-30X y MJV y MJV + 10i? DX + AR MN D m D n X + 12D M RD M X - 30X 2 i? 
+ \2X OR — hX R 2 + 2X R MN R AIN — 2X R ABMN R AB ^ IN 

Here I is the N x N identity matrix and Y MN is the matrix- valued quantity defined above in 
terms of the commutator to two covariant derivatives. 

Bulk Counterterms for Spheres 

Since our applications are to compactifications on spaces which are spheres, it suffices to 
specialize the general results of the appendix to these simpler background field configurations. 

Consider therefore 6D spacetime geometries which are the product of 4D Minkowski space 
with a maximally-symmetric 2D manifold: 

ds 2 = v, dzMx" + g mn dy m dy n , (A.9) 

where maximal symmetry for the 2D metric implies R mnvq = \R{gmp9nq - gmq9np), Rmn = 
\ Rgmn and D m R = 0, and so R mn p q R mnpQ = 2R mn R mn = R 2 . We also take any background 
scalars to be constants, and allow only a single background gauge field to be nonzero, and 
take it to be proportional to the 2D volume form: F mn = f e mn , for some scalar /. 



conventions [2£( ], which for the Riemann tensor agree with those of ref. |26[, but disagree with this reference 
by a sign for the Ricci tensor and scalar. 
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With these choices the only nonzero components of 'Y'mn lie? in the 2 dimensions, Y rnn ^ 
and all of the curvatures are covariantly constant. The coefficients a\ through 03 simplify 
considerably, reducing to |2£ 



a = I (A.10) 
01 = --R-X (A. 11) 

a 2 = — R 2 + -RX + -X 2 + — Y mn Y mn (A.12) 
^ 60 6 2 12 y ' 

11 1 1 

„ _ _ p3 ym -yn V 1 - __ f? V V mn Y V vmn 

— 630 30 n m 40 mn 12 mn 

--X 3 - —X 2 R- —X R 2 . (A.13) 
6 12 60 v ' 

Here / is the N' x N' identity matrix, and X and Y mn are the N' x N' matrix-valued quantities 
defined above. These expressions may be used to compare with the bulk part of the ultraviolet 
divergences encountered in the explicit calculations of the main text. 

Scalars, spinors and gauge fields 

This section collects the results for X and Y MN , and so also for the ultraviolet-divergent parts 
of the one-loop action, for several kinds of particles in 6 dimensions. Attention is restricted to 
those fields that do not mix appreciably with the gravity sector, and this means in particular 
that any gauge fields considered cannot be those whose background flux stabilizes the extra 
dimensions. 

Scalars 

For spinless fields we begin with the following general scalar-field action involving two deriva- 
tives or less, 



S = - d 6 x^ 



X - g MN Gij D M &D N & + V + ^UR+^H F a MN F^ IN 



(A.14) 



for a collection of iV real scalar fields, coupled to a background metric, g M N, and gauge 
fields, A a M . The functions U, V, W and the target-space metric, Gij, are imagined to be 
known functions of the $*. The background-covariant derivative appropriate to this case is: 

D M & = 8 M & - iA a M {ta) 1 ^ , (A.15) 

where the matrices (t a ) l j represent the gauge group on the scalars. 

To compute the one-loop quantum effects of scalar fluctuations we linearize this action 
about a particular background configuration, ip l , according to: $ l = (p l + , where d M ip = 0. 
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Expanding the classical action to quadratic order in <p l reveals the kinetic operator A l j, which 
is given by 

A^-^D + A 1 ,, (A.16) 

with X l j given by 

V hj {<p) + l -RU k] (tp) + \ F a MN F™ N U kj {ip)] . (A.17) 



X i. = G ik 



In this last expression the subscripts on U, V and % denote differentiation with respect to the 
background field ip l . Specializing to the simple 2D geometries and Maxwell fields discussed 
earlier, these simplify to X l j = G lk [Vkj + \RUkj + \f 2 1-Lkj\ and Y mn = -igfQe mn , where 
g is the gauge coupling and gQ = t a is the hermitian, antisymmetric charge matrix for the 
background gauge field. Notice that these imply Y mn Y mn = -2g 2 f 2 Q 2 and Y m n Y n \ Y l rn = 
0. 

Since scalars transform trivially under Lorentz transformations, d = and so N' = N. 
The ak are then N x N matrices, with the trace of ao through a 3 given by 

TV 

tra = N, trai = -— R-trX, (A. 18) 

and 

tra 2 = ^R 2 + ^RtvX+^tvX 2 -^~g 2 f 2 tTQ 2 (A.19) 

tr ° 3 = ~^k RZ + k R ~ 92f2iTQ2 + \~ 92f2iT{XQ2) 

--trX 3 - — RtiX 2 - — R 2 trX. (A.20) 
6 12 60 y J 

These give explicit functions of if once the above expression for X is used. 
Fermions 

For 6D massless spinors, ip a with a = 1, N, we take the following action 

S = -J dW^G^M^W, (A.21) 

where Tf) = e A M j A D M with 

D M r = d M r - utflAB^ - iA a M t a ^ , (A.22) 

with j A being the 6D Dirac matrices and e A M the inverse sechsbein, ^j AB = 5(7^, 7s], and t a 
denotes the gauge-group generator acting on the spinor fields. Since 6D Weyl spinors have 4 
complex components their representation of the 6D Lorentz group has d = 8 real dimensions. 
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The differential operator which governs the one-loop contributions is in this case If) = 
€a m ^ a D m an d go m order to use the general results of the previous section we write (assuming 
there are no gauge or Lorentz anomalies) logdet^) = \ log det(— fl 2 ), which implies 



i£ 1/2 = ^Trlog|> = ^Tr log (-^ 



= ^Tr log (-D -\R+ \l AB F a AB t^j . (A.23) 

This allows us to adopt the previous results for the ultraviolet divergences, provided we divide 
the result by an overall factor of 2 (and so effectively d = 4 instead of 8) , and use 

X = -\R + \l AB F a AB t a . (A.24) 



Similarly, we find 
and so 18 



Ymn — 2 Rmnab^I iF Mnta i (A. 25) 



MAT a bMN ^ ABM N 

Tr i/2[YmnY mn ] = —Atv l i 2 {t a tb) F^ IN F MN — — R ABMN R 



N 



= -8g 2 f 2 tr 1/2 (Q 2 )--R 2 , (A.26) 

where the second line specializes to the simple spherical 2D geometry and background gauge 
fields discussed earlier. 

Keeping explicit the sign due to statistics, and dropping terms which vanish when traced, 
this leads to the following expressions for the divergent contributions of ./V 6D Weyl fermions: 

(-) F Ti 1/2 [a ] = -AN, (-) F Tr 1/2 [oi] = ~R (A.27) 

(-) F Tr 1/2 [a 2 ] = ^R 2 - | g 2 f 2 tr 1/2 (Q 2 ) (A.28) 

(~) F Tr 1/2 [a 3 ] = -JL& + l?fRtr 1/2 (Q 2 ) . (A.29) 

Gauge bosons 

For N gauge bosons, A a M , with field strength Tl {N and a = 1, N, we use the usual Yang- 
Mills action 

S = -J d 6 *,/^ \ H(<p) T a MN K N , (A.30) 



18 We adopt the convention of using Tr[...] to denote a trace which includes the Lorentz and/or spacetime 
indices, while reserving tr [...] for those which run only over the 'flavor' indices which count the fields of a given 
spin. 



-47- 



expanded to quadratic order about the background fields: A a M = A a M + 8A a M . For an appro- 
priate choice of gauge the differential operator which governs the loop contributions becomes 

A aM bN = -S a b 5 M N D + X aM bN , (A.31) 

with 

X aM bN = -R M N 5 a b + 2i{t c )\F CM N , (A.32) 

where t c here denotes a gauge generator in the adjoint representation. 

The dimension of the 6- vector representation of the Lorentz group is in this case d = 6. 
We therefore find Tr v I = 6N, Tr v (X) = -N R, and 

Tr V (X 2 ) = N R MN R MN + 4 C(A) F« N F* IN (A.33) 
Tr V (Y MN Y MN ) = -N R ABMN R ABMN - 6 C(A) F a MN F^ , (A.34) 

where C(A) is the Dynkin index for TV" fields in the adjoint representation, defined by 
tr (t a t b ) = C(A)5 ab . The subscript l V J in these expressions is meant to emphasize that 
the trace has been taken over a vector field (as opposed to the physical spin-1 field, including 
ghosts). 

These expressions suffice to compute Tr v [afc], for the vector field. Once specialized to 
the spherical geometries and background gauge field of interest we find 

(-) F Tr v [a } = 6N, (-) F Tr v [a x ] = (A.35) 
(-) F Tr,[a 2 ] = ^fl 2 + 35 2 / 2 tr(Q 2 ) (A.36) 

(-f Tr v [a 3 ] = ~ R 3 + \ g 2 f 2 Rtr (Q 2 ) . (A.37) 
1U5 5 

To this we must add the ghost contribution, which consists of iV complex scalar fields 
having fermionic statistics and transforming in the adjoint representation of the gauge group. 
The contributions to the may be read off from our previously-quoted expressions for scalar 
fields in the special case X = 0. For such fields we have 

(-) F Tr gh [ao] = -2N , (-f Tr gh [ ai ] = y R (A.38) 
(~) F Tr gh [a 2 ] = - ^ R 2 + \ g 2 f 2 tr (Q 2 ) (A.39) 
(-) F Tr gh [a 3 ] = ^ R 3 - 1 g 2 f 2 Rtr (Q 2 ) . (A.40) 

Summing the contributions of eqs. ( A.35| ) and ( |A38 ) gives the contribution of physical 
6D massless gauge bosons: 

(-) F Tn[a ] = AN, (-) ir Tri[ai] = jR (A.41) 
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(-) F Tr 1 [a 2 ] = ^R 2 + f~g 2 ftr 1 (Q 2 ) (A.42) 
2/V 7 

(-) F Tr 1 [a 3 ) = -—& + - ~g 2 f 2 Rtr 1 (Q 2 ) . (A.43) 
B. Sums and zeta functions 

This appendix has two purposes. The first section provides a more secure theoretical foun- 
dation for many of the zeta- function calculations of the main text; while the second section 
provides formulae for several useful sums encountered in the calculations. 

Justifying zeta function regularization 

The purpose of this section is to compute the sum 



where 



d^(u,t,a ) :=^^g(27r&), (B.l) 

j=0 fc^O 



/oo 
dx exp[-i(j + w\x\ + a ) 2 ] e~ iqx , (B.2) 
-oo 



without resorting to the use of zeta function regularization. Pursuing this exercise will show 
that the zeta function approach used later is valid. 
Computing this Fourier transform gives 



2oj V t 

where 



(B.3) 



j := Vt{j + oo) , q := ^= ■ (B.4) 
Since G(—q) = G{q), we can write 

oo oo 

d rcs (oo,t,a ) = 2^^g(2mfc). (B.5) 

j=0 k=l 

When evaluating G{q) for non-zero q, it is important to note that the overall factor of 
exp(— q 2 / (4oj 2 t)) in eq. ( |B.3| ) will exponentially suppress any quantity multiplying it that does 
not diverge in the t — > limit. Since we are only interested in the t — > limit on this sum, 
we can therefore approximate G(q) as follows: 

G(q) ~ --/e-« 2 [e 2i ^erf(j + iq) + c.c. 



2oj V t 
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Uly/i 
1 



dxe x + 



j+iq 



dxe x + c.c. 



a; 



ie 2 ^~ / due-^ u+3)2 + 



c.c. 



(B.6) 



(Here, we have taken u = — i(x — j).) It is also helpful to make the following substitutions: 



Oil) 



v 7 * 



due u 3 sin[2j(g - u)] 



dy e 



g 2 (i-y) 2 



e 3 sin(2jqy) 



(B.7) 



where y := (q — u) jq. Since all of the j-dependence is in the bracketed term within the above 
integral, we can consider first carrying out the j-sum and then performing the integral over 

y- 



y<g(2irk) = -e-~ k2 / t 

J Cdt 

3=0 



dy e 



k 2 (l-y) 2 /t 



J2e~ ]2 sm(2fky/Vt)) (B.8) 



d=o 



:=R(ky,t) 

where we have introduced the notation k = irk/co. Calculating R(ky,t) can be done as before, 
using the Euler-Maclaurin formula in eq. ( |3.45| ), Defining 



g(x + a ; ky, t) := e - t{x+ao)2 sin (2(x + a )ky) 



this gives 



R(ky,t) 



D; 



dxg(x + a ;ky,t) - ^2~^-g^ l) {a ;ky,t) 



(B.9) 



(B.10) 



i=l 



where i/*" 1 ** (ao; ky, t) denotes the (i— l)-th derivative of g(x; ky, t) with respect to x evaluated 
at x = ao- The integral can be simplified as follows: 

-ao 



dx g(x + ao; ky, t) 



dx g(x + ao;ky,t) + / dx g(x + ao;ky,t) 



"0 



dx g(x;ky,t) — / dxg(x;ky,t) 



dx g(x; ky/Vt, 1) 



where 



D(w) ■= dxe x sin(2xw;) 



(B.ll) 



(B.12) 



is known as Dawson's integral - a real, analytic function with known series expansions about 
w = and w — > oo. The second integral can be reliably Taylor expanded in its upper 
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limit, since the t-dependence of the integrand is contained within the sine function (which is 
guaranteed not to diverge). Furthermore, since the integrand is an odd function of x, only 
even powers of \/taQ will appear. The coefficients in such an expansion will be related to 
derivatives of g(x; ky,t) (with respect to x): 



a _ w 

dxg(x;ky/Vi,l) = 



i=l 



(2i)! 



(tai)\ 



(B.13) 



Given all of this, we see that the sum over Q(2irk) can be divided into three contributions: 



j=0 



where 



Vt 



G l (k,t) + G 2 (k,t) + G 3 (k,t) 



Gl (jb,t) = JL e -»»/« / dye ^-y)Vt D ih 



(B.14) 



G 2 (M) 



fc \- (tag)* a /t 

(20! 7o 



dye 



Vt 



k 2 {i-y?/t n {2i-i) ( r,. 



(B.15) 



g^(0iky/Vt,i) (B.16) 



Si 



8=1 



(B.17) 



The small-i limit of these expressions can be obtained by making extensive use of the definition 
of the Dawson integral in eq. ( B.12| ), as well as its asymptotic expansion 



1 1 

D(w) ~ h 



15 



2w 4w s 8w 5 16w 7 



w 



(B.18) 



(This is most easily done using computing software.) And, of course, any term in the fi- 
nal result that is suppressed by a factor of exp(— k 2 /t) can be safely dropped in this limit. 
Performing this expansion yields 

' 2 3t 3 



Gi(M) 



t t z 
W + W + 16fc 6 



+ ... 



tai 



t 3t z 



2„4 



3t 



(B.19) 



In the above, the . . . indicate terms that will contribute to Sq(u, t) at 0(t 3 ) or higher, and so 
will give no contribution to the Casimir energy. 

Finally, we can substitute this result into eq. ( jB.5[) and, with the use of 

r2 ^4 ^6 



El /UJ\n , . 7T j TT , . 71"° 

F =(-) «») with <(2) = -, (i g).— 

k=l 



£ 

k=l 



1 1 

W ~ 6 



^ A; 4 90 ^ ' t.6 Q4 



fc=i 



fe=l 



A: 6 945 



we find that 



k=l 

oo 

4^G 2 (M) 

k=l 

oo 

4^G 3 (M) 



ojH 2 ojH 3 



+ 



+ 



6 180 1260 



+ ... 



a1uj 2 t 2 ciqUjH 3 (XqUjH 3 



k=l 



6 
1 

~2 



60 12 

aou 2 t 2 aQUJ 2 t s i aQUjH c 



+ 



30 



w 2 i 2 aluoH 3 uH 3 
— + 



-(0) 



6 
1 

"30 



60 



u 2 t 3 

12 



and so 



oo / oo 



d re s(w,t,a ) = 2 x ^ Y, Q ^ k ) 

k=l \j=0 



4 



oo 

[G 1 {k,t) + G 2 (k,t) + G 3 (k,t) 



k=l 










1 




[(■ 








~36 




+ 


[(- 



1 1 2\ 2 1 4 

+ gO0 - gO J U + — W 



1 



180 

1 



1 2 1 q ] \ ' ' 

H a n an H «n w 

360 12 6 12 



1 1 



1 



2 I , ,4 



H on a-n ^ 

360 60 60 ' 



1 



1260 



-a; 



as found later, using the more efficient zeta-function regularization method. 
Some useful sums 

Next, we will explicitly evaluate the small- 1 limit of two often-encountered sums. 



(B.22) 



(B.23) 
(B.24) 



(B.25) 



(B.26) 
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The sum c(t, ao) 

This sum is defined as follows: 

oo 

c(t, a ) := J2 e M~t(j + ao) 2 } . (B.27) 
j=o 

The Euler-Maclaurin formula states that, for any analytic function f(x), we can write 

J2fU)= dxf(x)-J2^f {i - 1] (V) (B.28) 
,=o ^ i=i *■ 

where f^ l ~ 1 \x) denotes the (i — l)-th derivative of /(x) and where the B{ are the Bernoulli 
numbers. (The first few are B x = -1/2, B 2 = 1/6, B 4 = -1/30, and B 2i+ i = for i > 1.) 
Identifying /(x) = exp[— t(x + ao) 2 ], the first few derivative terms are 

/(°)(0) = e- ta ° , /«(0) = -2ta e- ta ° , /( 2 )(0) = {-2t + 4t 2 a 2 ) e~ ta o 

(0) = (12t 2 a - 8t 3 af]) e~ ta o , (0) = (12t 2 - 48i 3 a 2 1 + 16t 4 ag) e~ ta o , (B.29) 

and /( 5 )(0) ~ 0(t 3 ). (We only need powers of t smaller than 3, for reasons discussed in the 
text.) The integral can be written in terms of the error function, and has the following series 
expansion: 



j™dxe-^ + ^ 2 = ly^(l-erf(a Vt)) 



Upon Taylor expanding the derivative terms, we find that 

c(t, ao) = — ^p- + c + ci t + c 2 t 2 + 0(i 3 ) (B.31) 

V* 

where 

-1/7F 1 1 1 2 1 o 

c -i/2 = ~2~> c := - - a , ci := -a - -a + -a , 

c 2 := ^a - ^ag + ^a 4 - ^a[> . (B.32) 
We will usually need a "tilded" version of the above sum. This is defined in the following 

way: 

c(t, a , r) := e Tt c(t, a ) = + c + c 1/2 Vt + ci t + 5 3/2 t 3/2 + c 2 t 2 + 0(t 5/2 ) . (B.33) 

In terms of components, 

C-1/2 = c_i/2 , c = c , c 1 / 2 =c_ 1/2 r, ci = ci+c r, 

c 3/2 = ^r 2 , c 2 = c 2 + Cl r + |r 2 . (B.34) 
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The double sum d(uo,t,ao) 
This sum is defined as follows: 

oo oo 

d(oj,t, ao) := exp[— t(j + uj\n\ + ao) 2 ] . (B.35) 

j=0 n=— oo 

The j'-sum can be performed as before, but the n-sum will pose a difficulty in the uj,t — > 
limit, where it converges very slowly This is because the sum in fact diverges in this 
limit. However, by using Poisson resummation, we can compute these divergent contributions 
efficiently. 

With this in mind, let's write 

d(u, t, a ) = d sing (uj, t, a ) + d rcg (oj, t, a ) (B.36) 

where we require that d rcg (uj, t, ao) contain all terms that vanish as u,J — > 0. Since the 
calculation of d sins (u, t, ao) is more subtle (for the above-mentioned reasons), we will compute 
it first and then apply a more cavalier approach to d reg (w, t, ao). 

Poisson resummation and d smg (LO,t, ao) 

A method for finding the behaviour of d smg (oj, t, ao) as cu, t — > is to use Poisson resummation. 
Defining g(n) := exp[— u> 2 £(|n| + bj) 2 ] where bj := (j + ao)/u, we can exchange 

oo oo 

E 9(n) = E (B.37) 

n=-oo fc=-oo 

where 

/oo 
dxg(x)e-^ x (B.38) 
-oo 

is the Fourier transform of g(x). For example, in the special case where bj = (i.e. j = ao = 
0), we find that g(n) = exp(— oj 2 tn 2 ) and so 



since 




(B.39) 



(B.40) 



This special case is sufficient to illustrate the value of Poisson resummation: a sum that 
converges slowly in the to, t — > limit is transformed into one that converges quickly in this 
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same limit. In fact, the part of the n-sum that diverges in this limit gets mapped entirely 
onto the k = term in the /c-sum. Therefore, we can write 



d si ^(u,t,a ) =J2 0(0) (B.41) 
j=0 



where now, for bj ^ 0, 



0(0) = f°° dx e-^'A+bj? = - erf (uVibj)] . (B.42) 



This sum can easily be evaluated in the same way as in section [B| using the Euler-Maclaurin 
formula. Integrating and differentiating Q(0) as prescribed gives 

1 / — 00 

C g (u,,i,a ) = -VlE[ 1 - erf (^' + a o))] (B.43) 
= °=±- + ^2J1 + d s mg + < ng t + cf 2 ms t 2 + 0(t 3 ) 

t yt 



in the small-i limit, where 

^f:=l, d*^ := ^ Q - oo) , df« := 1 Q - ao + og) , (B.44) 

,sing ._ j_ /_1 «| _|_ ^0 _ ^oA ^sing ._ 1_ f _}_ a , a a , a 

1 ' u V180 6 3 6 J ' 2 ' w V1260 60 12 10 30 
As promised, each of these contributions diverges in the limit where u,t — > 0. 

Calculating d TCS (uj, t, ao) using a zeta function approach 

First, we should mention that the use of zeta function method here is not necessary. In fact, 
we could continue as before, and compute d reg (uj,t,ao) using Poisson resummation: 

oo oo 

d ieg (oj,t,a ) = 2^J] G(2irk). (B.45) 

j=0 k=l 

Here, we have identified d reg (u, t, ao) as the left-over k ^ terms in the Poisson sum, 



eq. ( B.37 ). (The factor of two arises because Q(q) is even: g(—x) = g(x) =>• Q(—q) = Q(q).) 
This more "honest" approach is taken in Appendix [B[ However, since this calculation is 
significantly less tiresome if we take advantage of zeta function regularization, we will present 
this version of the derivation herein. (The two calculations give the same result, as can be 
seen by comparison with the result in Appendix [B].) Much of the simplification arises because 
we can take advantage of the previous result for c(t,ao). 
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Since we are guaranteed to have captured all of the divergent terms (as u>, t — > 0) in 
d smg (w, t, oo), computing the full sum in a naive way — one that does not capture these 
divergent terms — will give us d reg (cj, t, <Xq). The zeta function method is just such a naive 
approach! Let's start by interchanging the j- and n-sums in eq. ( jBjgp and write 

oo oo 

d reg (co, t, ao) = ^2 c(t, w|n| + ao) = c(t, ao) + 2 c(t, wn + ao) . (B.46) 

n=— oo n=l 

Since, the zeta function is defined as 

oo 1 

CW = E - > (B.47) 
we have that Y^=\ 1 = Ylri=i l/( n °) = C(0) = — 1/2 and so we can write 

oo 

ef eg (w, i, oo) = 2 ^ (c(t, wn + a ) - c{t, a )) . (B.48) 



n=l 



Given the result in eqs. ( B.31 ) and ( B.32 ), we find that the summand — expanded out in 
powers of t — is given by 

Ac := c(t, um + a ) - c(t, a ) = ^ f Ac r (B.49) 

r 

where r G { — 1/2,0, 1,2} and where 

/l n\ ( 1 \ n n 1 H 1 

Ac_!/ 2 = , Ac = -urn , Ac\ = I - - a + a j uju + - + a j ui n + 3^ n 

**=(^-4 + 4-4)»n+(-% + ^- 4) - 2 n 2 (B.50) 



m 2 2 y v 2 2 

+ (-1 + a - ag) a,V + Q - f ) ,V - lu, 5 n 5 . 

Each term in Ac proportional to n s will get a corresponding factor of C(~ s ) once the sum 
over n is performed. (Here, £(-1) = -1/12, £(-3) = 1/120, £(-5) = -1/252 and C(-2s) = 
for s > 1.) The end result is 

d reg (w, i, a ) = d v eg + i + d r 2 eg t 2 + 0(t 3 ) (B.51) 

where 

6 ' 1 V 36 6 6 / 180 ' y J 

,res ( 1 «0 a O a o\ /l a a nA q ^ 5 



2 V 360 12 6 12 / V 360 60 60/ 1260 



- 56 - 



As promised, the above result vanishes in the uj, t — > limit. 

Combining the results for d sme (oj, t, ao) and d reg (w, t, ao), we find that 

d(u, t, ao ) = ^ + +d + d 1 t + d 2 t 2 + 0(t 3 ) 



(B.53) 



where 



d-i 
d, 

d 2 



i , fi 



J 1 ^ 2 ^ 

d o = - [ ~ ~ a + a + — 
co \6 6 

4 1 



1 Q o , a o Q o , f_ J_ . a o a oV.2 , ^ 
180 6 3 6 \ 36 6 6 / 180 



(B.54) 



n 2 n 4 n 5 n 6 



1260 60 12 10 30 

I __L + ^°_ , 
1 360 60 60 ' " 



1 , a o a o , a o , ... 
360 12 6 12/ " 



CO 



1260 



As with c(t, ao), we will often need a "tilded" version of the above sum. This is defined 
in a similar way as before: 



d(uj,t, a,Q, t) := e T d(uj,t,ao) 



(B.55) 



= d_i + d_^2 + ^ + ^ + ^ t + t s/2 + J 2 t 2 + ^5/ 2) _ 



3/2 1 



In terms of components, 



d_! = d_i, d_ 1 /2 = d_ 1 / 2 , d = d + d_ir, d 1 / 2 = d_ 1 / 2 T, (B.56) 



Ji = di + d r+ ^— !-r 2 , d 3/2 



-1/2 2 J t I i ,^02, 3 

— A_ r ^ 5 d 2 = d 2 + di r + — r z + — r° 



C. Spectra and mode sums 

This appendix computes the KK spectra and mode sums that are quoted in the main text, 
for spins zero, half and one. 

Spectrum and Mode Sum for the Scalar Field 



This appendix computes the spectrum A^ n and the corresponding small-t limit of the mode 
sum 

S.(t):=£ e -**3» = E^ + "- (C- 1 ) 

j,n i=— 1 

for a 6D charged scalar (ft, with mass m and charge qg, on the rugby ball. 
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Eigenvalues 

The scalar equation of motion is 



! D M D, 



m 2 ^<p = 



where D M is the covariant derivative, and where 



ds 2 = ri^dx" dx v + r 2 (d0 2 + a 2 sin 2 9 dp 2 ) 
(N - 



qA m dx m = qA^dip 



■(cos 6- b) + b$ t 



dip . 



(C.2) 

(C.3) 
(C.4) 



In the above expression for A^, the variable b is used to distinguish two patches of the gauge 
potential: b = +1 (—1) corresponds to the patch encompassing cos6> = +1 (—1). To verify 
that this gauge potential is a correct one, it is sufficient to check that: 



dA v 
d6 



fee 



M 

2qr' 

A <p (cosO = b) = b^. 

Z7T 



ar sm ( 



(C.5) 
(C.6) 



(Recall the definitions N = uj(N - $), & b = qA h /(2ir) and $ = from the text 

Expanding out this equation of motion using the ansatz (f>(x 



e ik 4 -x ®(Qyn V gives 



d 2 e 



d9< 



+ cot I 



de 

~d9 



+ 



(n-qA^f 
a 2 sin 2 9 



e = o 



(C.7) 



where Xj n = —r 2 (kf + m 2 ) are the required eigenvalues (in units of the KK scale 1/r). 
This is just the result for the case without a background field, but with n y in q Aip . 
This difference is expected given that the wave equation depends on the azimuthal covariant 
derivative D v = d v — iqA v = i(n — qA v ) (since our ansatz is that <f>(x) oc e 



inip\ 



The standard approach from here is to change variables to x := cos 9, giving 



(1 



d 2 @ 

dx 2 



de 

2x— + 

dx 



[n 



e = o, 



(C.f 



and to let 6(x) 



a 2 (l — x 2 

(1 — x 2 ) y f(x) for some convenient power y that removes the singular 
behaviour of ( C.8 ) near x = ±1 9 = 0, 7r). However, since the gauge field is defined 
differently at each pole, it will help to treat the divergences at each pole separately. Therefore, 
we will instead take 

e(x) = (i-x)y(i + xyf(x), 

which yields 



d 2 / 



df 



(1 - x 2 )-^ - 2[(1 + y + z)x + y - z]-± + 



dx 



Xj r 



u 2 (n - q A^) 2 + k(x, y, z) 



(1 



x' 



e = o (c.9) 
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where to = a 1 and 

k(x, y,z) := y + z - (y - z) 2 - 2(y 2 - z 2 )x - (y + z)(y + z + l)x 2 . (C.10) 

Requiring that the numerator [(n — gA/p) 2 + fc] be proportional to (1 — x 2 ) gives the conditions 

y = -\n + {N- $)e_ 6 - b$ b \ , z = -\n-{N - $)e fc - 6* 6 | (Cll) 

where e + = 1 and e_ = 0, and the resulting regular ODE is 

N 2 



{1 ~ x2) ^~ 2[{1 + y + z)x + y ~ z] ^ + 



X jn + — -(y + z)(y + z + l) 



f(x) = 0. (C.12) 



In passing, it should be mentioned that it is equally valid to take y and/or z to be definitely 
negative, rather than positive, in eliminating the singularities from eq. (|C.8| ), However, in the 
absence of couplings of the scalar to the brane, we should discard these solutions since they 
cause (j>(x) to diverge at one or both of the branes (in conflict with the boundary conditions). 

It is instructive to consider the special case in which q = and u = 1. In this case, we 
find that y = z = \n\/2 (thereby justifying the usual choice of Q(x) = (1 — x 2 )' n l/ 2 /(x)) and 



1 " i 2 )7T " 2 ( X + M^T" + IV " M(M + !)] /(*) = • 



( _1_ 2(1 + W)I £, .„„ 



This is the standard differential equation that is satisfied by the |n|th derivative of a Legendre 
polynomial (so long as \j n = (j + \n\)(j + \n\ + 1) for some j > 0, allowing the corresponding 
power series solution to terminate). Therefore, 0(x) can be any (linear combination) of the 
associated Legendre polynomials and <j)(x) are (linear combinations of) spherical harmonics. 
It is conventional to define £ := j + \n\ so that A = tit + 1) with degeneracy (2t + 1), 
representing the number of different ways of choosing combinations (j, n) that give the same 
t. 



In the general case, the condition for termination of the power series solution to eq. ( C.12 ) 

is 

M 2 

A jn = (j + y + z)(J + y + z + 1) - — . (C.13) 

Mode Sum 

From here, we wish to compute the small-i coefficients of S s (u, N, t), where 

oo oo 

s s (u, n, $ 6 , t) = e m2 ' A E ex p [ - ^ + y + + y + z + ^ ( c - 14 ) 

j=0 n=— oo 

oo oo 

= eTt E E ex P i-t(j + y + z + 1/2) 2 ] . (C.15) 

j'=0 n=— oo 
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In the last line, we have introduced r := (1 + A/" 2 )/4. 

In order to use our previously-derived results for the sums 

oo 

c(t,a ,T) := e r *^exp[-t(j + a ) 2 } 

3=0 

oo oo 

d(w,t,ao,r) := e rt ^ exp[-t(j + w|n| + a ) 2 ] , 

j=0 n=— oo 



(C.16) 
(C.17) 



we will first have to do some massaging. We will assume throughout that |<3?fc| < 1, since any 
brane-localized flux larger than this can be absorbed (one integer at a time) into the bulk 
flux, N. Also, we will choose to work with the north (cos# = +1) patch of the potential, for 
which 



y + z = — |n-$ + | + — |n-7V + <I>- 



(C.18) 



To see the result for the south (cos# = — 1) patch, where 



y + z = - \n + N- $ + | + -|n + $_| , 



(C.19) 



we can simply make the replacements N — > —N and $6 — > — in the final result (since the 
spectrum only depends on the combined quantity y + z). Of course, such a transformation 
should be a symmetry of the result, since it is just a gauge transformation. 

It is helpful to first consider the case when N > and break up the n-sum in eq. (C.15) 
into three parts: 



S e (u, N, $ 6 , t) = Si + \u, N, $ 6 , t) + N, * 6> *) + S(-\u, N, $ 6) i) 



(C.20) 



where: 



for <J>_ > 0, 3>+ > 0, 



oo oo 

j=0 n=N 
oo AT-1 

Ss (0) =e rt ££ ex P 

j=0 n=l 

oo — oo 

j=0 n=0 



. / | j + |( n _$ + ) + |( n _ AT + + ! 



./ j + |( n _ $ + ) + |(at - - n) + ~ 



-t ( j + -($+ - n) + -(JV _$__„) + - 



- 60 - 



• for > 0, $ + < 0, 



oo oo 

j=0 n =N 
oo AT-1 

j=0 n=0 

oo — oo 

s t ] = ^E E «p 

j=0 71= -1 



-/ ( i + |(n - <D + ) + |(AT _ _ n ) + 1 



woj 1 
-t ( j + -($+ - n) + -(iV - - n) + - 



for $ < 0, $+ > 0, 



oo oo 

^W'E E -p 

j=0 n=N+l 

oo N 

5 s (0) = e rt EE ex p 

j=0 71=1 

oo — oo 

s t ] = eTt E E ex p 

j=0 n=0 



_/ ( i + |( n _$ + ) + |( n _iv + $_) + I 



/ ( ./ + |(n - <&+) + |(iV - d>_ - n) + i 



a; , oj , 1 

M i + g (*+ - n ) + 2 (iV " $ ~ " n) + 2 



• for < 0, $ + < 0, 



oo oo 

^W'E E «p 

j=0 rt=JV+l 

oo AT 

^ (0) = e rt EE ex p 

j=0 n=0 

oo — oo 

^ } = ^EE ex p 

j=0n=-l 



, OJ , OJ , . 1 

M i + 2 ( n - $ +) + 2 (n " ^ + ^ + 2 



/ ( j + |(n - <&+) + |(iV - d>_ - n) + 1 



. / ( j + |($ + _ n ) + |(JV - <D_ - n) + 1 



If we allow ourselves to extend the notation e± to 

1 + sgn(x) _ J 1 , sgn(x) := x/|x| = +1 
2 | , sgn(x) := x/\x\ = —1 



(C.21) 



(recall that, before, e& = (1 + i>)/2 = 1 or 0), these four cases (for which N > 0) can be 
succintly written as 



oo oo 



S ( s +) =e^ E ex P 

j=0 7i=e_$_ 



l+wiV\ 2 ' 
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= (N- 1 + x e T * E exp 

j=0 



, . u§> 1+ujN 



oo — oo 

r '"E £ exp 

j=0 n=-e_$ 



wA$ 1+loN 
-t [ j — con -\ 1 



(C.22) 



where A<E> := <E> + — <E>_. 

Performing a similar calculation for the case where N < 0, we find something slightly 
different: 

2" 



oo oo 



J=0 n=e$. 



wA$ 1 - wiV 
-t ( j + wn — + 



5(0) = { _ N _ 

j=0 



cj<1> 1 - wiV 



oo — oo 



St 



e-E E ex P 

jr =0 n=— e<5_ 



wA$ 1 - wiV 
■t j-wnH 1 



(C.23) 



Writing r] : = sgn(iV) = 7V/|iV|, we find that these expressions - for arbitrary N - take 
the form 

*s (+) = e-E E 



wA$ 1 + oj\N\ 
-tij+oon — + -!— ^ 



3=0 n =£-v<s>- v 

oo 

SP = (\N\ - 1 + e- v q>_ + e_„* + ) x e rt E ex P 

j=0 



-mj-*-2- + — 2 



oo — oo 



s s ( -W'E E -p 

j=0 n=-e_^ 4>1) 

Using the final rearrangement 



wA$ l+wlATl 
■i h- wn H 1 



(C.24) 



oo 1 / oo \ 

E/w = 2 E /(H)±/(o) 

n=e:p \n=— oo / 



(C.25) 



and writing a N := (1 + u\N\)/2, we find that the massage is complete: 

S s +) = 2 d(w,t,-a;A$/2 + a JV ,r) + ??sgn($_ T; )g(t,-wA$/2 + a JV ,r) 
S B (°) = [|7V|-r ? (sgn($ + ) + sgn($_))/2] c(t, -r?a;$/2 + 0jv , r) 
= \ \d(u, t, wA$/2 + a N , t) + r) sgn($„) c(t, wA$/2 + ffljv ,r) 



(C.26) 



This result is, by construction, invariant under N — > —N, &b — > — <!>_&. 
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Evaluating these expressions, we find, using throughout the definitions 
F b := |$ b |(l-|$ 6 |) , F^:=J2 F b, fW-.= F, G{x) := (1 - x)(l - 2x) (C.27) 

b 

that s s _ 1 = 1/uj, 



s s (w,iV,$ 6 ) 
sl(oj,N,<f> b ) 

s s 2 (uj,N, $ 6 ) 



1 

a; 

1 

a; 

1 

UJ 



180 24 18 v 7 



12 



^a> 6 G(|^|) + ^(l-15F( 2 )) 



(C.28) 



(£.29) 



11AA 2 / 1 AT 
+ 



504 720 V 90 144 / 



+ 



i?(2) J?(3) ' 



1260 120 60 



(C.30) 



Spectrum and Mode Sum for the Spin-1/2 Field 

This appendix computes the spectrum A* n and the corresponding small-i limit of the mode 



sum 



-tA' 



(C.31) 



£ ii n i=— 1 

for a 6D charged spin-1/2 field ip, with mass m and charge q, on the rugby ball. Here, the 
minus sign denotes fermionic degrees of freedom, and the additional sum is over spin states. 



Eigenvalues 

The spinor equation of motion is 



(J0 + m)V> = 



(C.32) 



where D M is the covariant derivative and where the metric and background gauge field, 
respectively, are 



ds 2 = r 2 (dfl 2 + q 2 sin 2 9 dp 2 ) + rj^dx^dx" 



qA M dx ' = qAp dip 



-(cos9 - b) + 6$ 6 



dip 



(C.33) 
(C.34) 



(as before). Here, q is the charge (in units of g) of the field ip under the gauge field. 
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We make use of the frame fields e A M and the spin connection Q A f to ensure that our 
Cartesian Clifford algebra (i.e. {r A ,r s } = 2i] AB ) conforms to the spherical geometry: 

JP^ = T A e A M (d M - l -J AB n A M B - iqA M \ V • (C.35) 

(Here, J AB = — i\T A ,Y B \/A are the Lorentz generators in the spinor representation.) More 
precisely, these are given by two separate patches (in the same way as the gauge potential): 



i 



m _ _ 



b sin<z> 

COS ifi Ha 

~ a sin 

cos ip 

sin 9 



V — 



I • CO 

bsimp — 



<5S , e a m = (so that e A M e B N g MN = rj AB ) , (C.36) 



K B = V AC (e c N r K MN e K B - e c N d M e N B ) (where e M A := {e A M y') ■ 



(C.37) 



(Recall that 6+1 (6 = —1) for the north (south) brane.) Given that T^ N = with the 
exception of = cot 9 and T e w = —a 2 sin 9 cos 9, we find that the spin connections mostly 
vanish as well, with the exception of SI 4 , 5 = a cos 9 — b = — $7^ 4 . 

In order to expand out the above Dirac equation, let's use the representation 



r 4 = 



o 7 M 

V7 M ) 

where 7^ are the usual 4D Dirac matrices: 



75 
75 



r 5 = 



7" = -i 



a" 
-a" 



75 



-i7°7 1 7 2 7 3 



-il 4 
il4 



I 2 
-I 2 



and ct m = (l2,<7j) = cr^ are the Pauli matrices. 

Substituting these definitions, we find that the covariant derivative becomes 



= + T 4 



+r 5 



cos ip „ 6 sin v> . . , . cos <p 



^ 5 



a r sin t 

cos 
a r sin ( 



(c^ - iqA^) + 



a r sin 9 2 
sin ip 



v a r sin 6> 2 



or 



IP 



7^ + P L 2 - P fl Oi 

x 7"^ + P i Oi-P«0 2 
where P L (i?) = (1± 7s)/2 and where 



Oi 



O2 = 



{-id v - qA v - Qf/2) 



a sin 6 



asm( 



(C.38) 



(C.39) 



(C.40) 
(C.41) 

(C.42) 
(C.43) 
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Since tp is a 6D Dirac spinor, we can decompose it into two 4D Dirac spinors or, equivalently, 
four 4D Weyl spinors: 

^=(M = (^ + M (C.44) 

\ ^2 / \lp2L+ 4>2R J 

where the 4D Weyl spinors satisfy 'y^ipiL = +ipi L and jsifiiR = —ipi R - In terms of these, the 
Dirac equation gives 



^d^j 2R + 2 ^2l + mipiL = , 'fd^iht ~ Oii>2 R + mifiiR = (C.45) 
-fdrfm + OiVu + W>2l = , t^Vil - 2 ^i B + mfc = (C.46) 



(where we have applied the projection matrix identities P 2 = P\ = 1, P l Pr = 0). From these 
expressions, it is clear that setting m = decouples ipi from ^2- In other words, for massless 
fermions there is a halving of the number of degrees of freedom. 
Decoupling (|045|) and flCiffi) gives 



(□ 4 + O2O1 - m 2 )^ = , (D 4 + Oi0 2 - m 2 )Vi« = 
(□4 + O1O2 - m 2 )V>2z, = , (D 4 + O2O1 - m 2 )^ = 



(C.47) 
(C.48) 



where D 4 = r}^ v d^d v is the 4D d'Alembertian. In what follows, we will quote results only for 
tpi, since a solution for ipi L (ipi R ) will also be a solution for Tp2 R (^21, )• To obtain the mode 
sum for a massive spin-1/2 field, we should remember to multiply the result for ip\ by 2. 
Using the ansatze 

^i L (R)=e ik '- x e L{R) (6)e m ^ (C.49) 

where n L and n H are integers satisfying n R — n L = b given eqs. ( Ejjp , (|046D (and because of 
the e* 1 ^ factors in Oi, O2)) we find that the spinor functions &l(r) must satisfy the following 
equations: 



d 2 e 



L +cot^ + 



d6 



uj~ n 



1/2 



i«i/2 



1 

4 



[( n l/2 - tfAp) cos ^ + QF0<p sin #] + 



sin 



Q L = (C.50) 



d 2 e fl de fl 



+ 



i«i/2 



where Fg v = dA^/dO (in our gauge, where Ag = 0), 

"1/2 : = n L + b/2 = n R - b/2 , 



(C.52) 



and where 2 := — (k 2 + m 2 )r 2 are the required eigenvalues (in units of the KK scale). 
These equations transform into one another when L ^ R and ni/2 —n\/2, A v o —A^. 
Since such a transformation is equivalent to a change of coordinates (in particular, the ex- 
change of north and south poles), we expect that the resulting eigenvalues will be the same 
when solving either equation, as needed. The exception to this logic is the situation where 
k 2 = —m 2 (i.e. when A f = 0); in this case, eqs. ( |C.45| ) and (U.46) partially decouple. In this 
case, we will find that only one of ipi L , ipiR wm have zero mode solutions, given some fixed 
value N ^ 0. 

Remarkably, after: 1) performing the coordinate substitution x := cos#; and 2) substi- 
tuting 

e a ( x ) = {i-xy°{i + x) z °u(x) (c.53) 

where a = +1 (a = — 1) denotes the left-handed (right-handed) field, we find that equations 
flOoOD and ( |OolD become 



(1 - x 2 )^-^- - 2[(y a + z a )x + y a - z a }-^- + 
dar dx 



1 aN 
4 + ^T 



(ni/ 2 - qA v ) 2 - aui(n\/2 - qA^)x + k(x, y a , z a ) + \ 



(C.54) 

fa = 



where we have borrowed the previous definition for k(x,y,z) in eq. ( |C.1C| ) . (We have substi- 
tuted for Fqm to obtain the term proportional to N in the equation above.) Requiring the 
numerator in the second line of eq. ( C.54j ) to be proportional to (1 — x 2 ) makes eq. ( |C.54[ ) a 
regular ODE and gives the following result for y a , z a : 



1 

Va= 2 



uj n 



H/2 + (N -$)e_b-&$ 6 ) -- 



a 



uj\ n 



l l/2 



(N-<S>)e b -b$ b )+- 



a 



(C.55) 



The resulting differential equation is 



(1-x 2 



dx 2 



2[(1 + y a + z a )x + y a - Za]^r~ 

ax 



+ 



j"l/2 



1 TV 2 

4 + — - {y<T + z*){y<r + z<r + ij 



(C.56) 



In passing, it should be mentioned that, if the goal is regularity of eq. ( |C.54 ), it is equally 
valid to take y a and/or z a to be definitely negative, rather than positive, in eliminating the 
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singularities from eqs. ( |O50| ) and ( |05l| ). However, we should discard these solutions since 
they cause for ip to diverge at one or both of the poles. 

In the case where the spinor does not couple to the background gauge field (q = or, 
equivalently, N = ^j, = 0) and the background geometry is that of a sphere (a = uj = 1), 
the above equations simplify as follows. Keeping in mind that \ni/2\ > 1/2, we have that 
Ul = Zr = |(K/2l - \) an d y R = z L = ^(|n 1/2 | + \) and so eq. (|C.54|) becomes 



(l-^-2[( 1 + | ni/2 |),-l]f + 



\f 

jnx/2 



I "1/2 1 + 



U = (sphere). (C.57) 



The power series solutions to these equations terminate as long as 

A k/ 2 = C? + K/2I + V2) 2 (sphere) (C.58) 

for some integer j > 0. If we define £ := j + |%/2| — 1/2, the above expression for A f simplifies 
to the familiar {t + l) 2 , as found in [21]. The degeneracy in this case is 2(£ + 1), which 
accounts for the number of different ways one can choose (j, n-11% ) to get the same I. These 
two formulations are equivalent. 

In the general case, the condition for termination of the power series solution to eqs. (C.56) 
(at the same order x 3 ) is that 

where j a £ {0, 1, . . .} (and j L ^ j R generally). 
Mode Sum 

The spin- 1/2 mode sum for ip\ is 

St (w, iV, t) = 2S t ^u, TV, $ 6 , t) + 2S fR ( W , iV, $ 6 , t) 

where 



IV M 2 

Ja + Va + Z a + - I - — 



(C.59) 



(C.60) 



S ttr (u>,N, $ b ,t) 



-2 x e^ 2 / 4 



00 00 

2 ex p 

j'=0 n=— oo 



-i(j + y<x + ^ + 1/2) 



(C.61) 



Here, a = L = +1 (a = R = —1) gives the result for the left-handed (right-handed) spinors. 
(Note that we have dropped the a subscript on the dummy indices j, n.) Recalling that 
n i/2 ■= n + ba/2, we can write y a and follows: 



Va 



uj 



UJ 



n + (N - $)e_ 6 - b $ 



<7 




. ~ 2 


~~ 2 



n-(JV-*)e6-6 (*6- j 



(7 s 







6J 




): 






n 




+ 2 


~~ 2 





n + (JV f(T - $ f<r )e_ h - h <l> 
n- (iV f(T -$ f(7 )e 6 -6$ 



fa 



fcr 



(C.62) 
(C.63) 
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where, in the second expressions, we have introduced 



N fer :=N-a, $f° 



f(l-I) (and%„:=E^ 



$-(7 1 



(C.64) 



(Later on, we will also use $g := (1 — oj 1 )/2.) Given these identifications, we can write the 
fermionic mode sums in terms of the previously-derived scalar result: 

S {(T (oj, N, d> 6 , t) = -2 x e^-d+^/^w, JV fff , <D f 6 CT , t) (C.65) 

where Af CT := w(iVf CT — 5>f CT ) = ftf — a. Expanding this out, we find (using the notation 
introduced in eq. ( |C.27| )) that = -2/oj, 



4 <T (^,iV f(7 ,^ CT ) 
4>,iV fCT ,<I> f 6 CT ) 



1 

a; 

1 

a; 



1 

6 ~ 

7 
720 



oN + 

_7V 2 
6 



+ AF fCT w 



(C.66) 



+ , 72 



1 (1 - 3ajV) 



6 



+ 



oj 3 (jV - a) 



6 



+ 1 — + Wfi 

V 720 6 L la J 



(C.67) 



31 



317V 2 



20160 1440 
oj 3 jV(1 - aAA) 



+ 



TV 2 



12 



2880 144 



;i-5aa 2 ) . 2 

-AF f(J w 



ri)]+ 



60 
7 



77V 2 



(1 - 3ajV + A/" 2 
12 



A[^] W + 



2880 1440 



60 



+ 



31 A[F f ( CT 2) ] A[F f ( CT 3) i 



+ 



+ 



20160 60 30 
where the Ffo b are calculated using and where we have taken 

A[X] := X -X 



^A[^G(|$n)(l + 3^)] 
b 

(C.68) 



(C.69) 

to ensure that A[X] -> as -> 0. For example, AF fe f<T = 1^^1(1 - |*^|) - (1 - w" 2 )/4. 
(Recall that oj > 1 for the rugby-ball.) It is important to remember that the derivation of 
these results has assumed < 1, which is distinct from the restriction |$f,| < 1 in the 

scalar case but of little difference near oj = 1. 

Since the 6D the chiral spinor ip\ is the sum of a left-handed and a right-handed 4D 
chiral spinor, its s\ coefficients are given by 



cre{±l} 



Jo 



(C.70) 
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(For massive 6D spinors, we must double this result.) Since we expect some cancelation when 
considering this sum of coefficients, let's work out the result. We must be careful, though, 
since there is some subtlety in this. For example, the spin-averaged quantity (AF fe f ) depends 
differently on for |<J>&| > <I>q := (1 — oj~ 1 )/2, as compared to when < $q: 



in 



-n 



I IrlS. I >> <T.f 



F b -%, \$ b \ > % 



(C.71) 



Both of these expressions have useful domains of validity: when considering small fluxes, <£;,, 
in the presence of some <£q / 0, it is the upper function which is of interest; in the case of the 
larger fluxes, or in the sphere case where oj = 1 and <l?g = 0, it is the lower function which is 
relevant. We should therefore treat these two situations differently. 
In the case where both |<J>&| < $q, we find that 



4(oj,N, $ 6 ) 
s[(oj,N,<f> b ) 



1 








M 




OJ 






1 


7 


ojjV® 


OJ 


360 ~ 


2 



OJ 



(C.72) 



+ 



36 6 E $b 



OJ 



OJ 



E^-^ + l^-sE^ 1 - 2 ** 



sI(oj,N, $ 6 ) 



1 

OJ 



31 ojN<5> 31 N 2 



360 
1440 



OJ 



b 

72 



(C.73) 



OJ 



w 



-OJ 



24 



V240 6 5 



+ 



31 



10080 



-E*i 



1 240 12 15 



(C.74) 



When both |$&| > $q, the small-t coefficients are given in terms of 



p b := sgn($ 6 ) 



$b 



$ b := w($ 6 - , p: 



F 6 := |$ 6 |(1-|$ 6 |), F:=Y,h etc. 



(C.75) 



as follows: 

OJ 



-- + 2F + 2oj 
3 



2oj 2 



(C.76) 
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s[(u),N, $ 6 



1 



1 F Af& pAf\ 2 lo a 
+ 9"¥ + — --6" w -45 



M 2 



(C.77) 



4(w,JV,$ 6 ) = - 



1 F® m P N Af ( & b F b ¥ b 



315 30 



15 



45 

1 /•' 



90 



6 



-AA 2 ( — - — ) I u; 2 + ( — — — — — + + — 
18 6 I 90 30 30 30 90 



315 



(C.78) 



(Note that we should take = when comparing these sets of equations when they overlap 
at |$»{,| = <1?q.) In the case of an uncharged, massive fermion on a sphere (i.e. when oj = 1, 
M = $>& = 0), these two sets of expressions both give 

8^(1,0,0) =-8, 4(1,0,0)=!, 4(1,0,0) = ^, 4(1,0,0) = J- (C.79) 

3 15 63 

in agreement with the result in 



Spectrum and Mode Sum for the Massless Spin-1 Field 

This appendix computes the spectrum A| n and the corresponding small-i limit of the mode 



sum 

j,n i=—l 

for a 6D massless spin-1 field A M on the rugby ball. 
Eigenvalues 

The spin-1 equation of motion is 

g MN D M F NP = 
where F MN = d M A N — d N A M , D M is the covariant derivative, 

{Dm) p a F NP = d M F^ P — T2 in Fq P — VfjpF^g — i(t ) p a A M F NP , 

and the metric is 

ds 2 = r 2 (d8 2 + a 2 sin 2 d dip 2 ) + r)^dx^dx v 

= r 2 (d8 2 + a 2 sin 2 6 dip 2 ) + 5 ij dx i dx j + 2dx"dz + . 



(C.80) 



(C.81) 



(C.82) 



(C.83) 



- 70 - 



In the last line above, we introduce light-cone coordinates x^ := (x 3 ± x°)/\^2 which are 
convenient for what follows. (Also, i, j run only over 1, 2.) More formally, we can introduce 
the (sub-) metrics g a p and g mn over the light-cone coordinates a,/3 G {+, — } and extra- 
dimensional coordinates m,n £ {9,f}, respectively: 

/oi\ (r 2 \ 

9af3 = , fan = , _ _ • C.84) 

\ 1 J \0 a 2 r 2 sin 2 6 J 

The full metric g MN is then composed as follows: 

9mn = 9mn © • (C.85) 

In what follows, we assume that the generator for the background gauge field A M is 
diagonal in the adjoint representation: (t A ) p a = q5°. Given this choice, it is convenient to 
suppress the group indices in what follows. We also choose the light-cone gauge, in which 

A- := As ~ A ° =0. (C.86) 
v2 



As before, the background gauge field is taken to be 

N- $ 



qA M dx =qA ip d(p 



(cos 6 - b) + b$ b 



(C.87) 



where N is an integer. Expanding out the equation of motion for P = — , we find that 

p = - : = - fidid-Aj - g~ + d 2 A+ - g mn D m d^A n 

= -d- {diA 1 + d-A+ + g mn D m A n ) . (C.88) 

This can be integrated to give the constraint 

F_ + = d _A + = -d^ - g mn D m A n . (C.89) 

This constraint tells us that A+ is completely determined (up to some initial conditions) once 
we have solved for Ak and A p . As we shall see, this constraint also simplifies the equations 
of motion for Ak and A p . 

Evaluating the equations of motion using the above constraint for P = k G {1,2} and 
P = p £ {6,(p}, respectively, we find that 

P = k : = d'Fik + g- + d-F +k + g + -d+F_ k + g mn D m F nk 

= (U A + g mn D m D n )A k (C.90) 

P = p : = d l F ip + g- + d^F +p + g+~ d + F_ p + g mn D m F np 

= (D 4 + g mn D m D n ) Ap + D p (g mn D m A n ) - g mn D m D p A n (CM) 
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where O4 := r}^ v d^d v = d l di+2 d-d+. Eqs. ( |C90 ) are exactly the same as the one encountered 
in the case of a scalar field, so the corresponding spectrum for A\ and A2 will be that of a 



scalar. Eqs. ( C.91 ), however, are coupled and our task is to decouple them. 

Given that the Christoffel symbols mostly vanish, with the exception of Tg^ = = cot 9 
and r 9 6 , = —a 2 sin 9 cos 9, we find that, using the ansatz 



Am(x, e, if) = e^- x ® m {ey n * , (c.92) 



Eqs. flC.91|) become 



(C.93) 
(C.94) 



d 2 e e 


+ cot 


Q d&g 


+ xe e - 


Kg 


d9 2 


9 de 


a 2 sin 2 9 


d 2 e^ 


— cot 6 




+ A9 V - 




d^ 2 


de 


a 2 sin 2 9 



where A := —k 2 r 2 are the required eigenvalues and where 



K e = [(n - qA v f + a 2 ] Og + 2»(n - qA 9 ) cot 9 9^ (C.95) 
K v = {n- qA^) 2 ®^ - 2i(n - qA v )a 2 cos 9 sin 9 ®e • (C.96) 

In order to algebraically decouple the K m terms, we should first perform a field redefinition 
such that the differential operators are the same for both equations. (Note that they are not 
currently, due to the sign difference in the cot 9dQ m /d9 terms.) 

With this goal in mind, we should (for convenience) switch to x := cos#, and substitute 

Q m (x) = (l-x)^(l + x)^f m (x), (C.97) 

where y m , z m are to be related in some specific way. Doing so gives 

(1 - X 2 )^ - 2 [(1 + Zg + yg)x + yg- Z g}^ + Xfg - = (C.98) 

dx z dx 1 — x L 

(1 - x 2 )^ - 2 [(z v + y v )x +y v - Ztp ]^ + X U - JS*- 2 = (C.99) 

for some K m , to be specified shortly. Equating the coefficients of df m /dx, we find that the 
required relations are 

yg = y v - 1/2 :=y, z e = z v - 1/2 := z . (C.100) 

Given these identifications, we find that the K m are given by 

Kg = (^(n-qA^) 2 + l + k{x,y,z)) f + 2iu 2 (n - qA v )x f v (C.101) 
K v = (u: 2 (n - qA^) 2 + 1 + k{x, y, z)) f v - 2i(n - qA^x fg (C.102) 
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where k(x, y, z) is defined as in eq. ( |C.1C| ). 

From here, we can see that the complex choice 



yields decoupled equations of motion for T (and F* ) : 



(1 -x 2 



d 2 F 

dx 2 



2 [(1 + y + z)x + y - z] — + f A - — ^ 



J" = 



(C.103) 



(C.104) 



where 



K T T := Kg + iuKy = [uj 2 (n - qA v f + 1 + k(x,y,z) + 2uj(n - q\)x] T . (C.105) 

Now that we have decoupled the equations of motion, we can eliminate the singular behaviour 
as x — > ±1 by requiring that Kjr be proportional to (1 — x 2 ). Doing so gives 



uj 

y = — 

y 2 



n + (N - $)e_ 6 - + - 



UJ 



yielding the spectrum 



j + y + z + - 



n - (N - $)e b - b$ b 

OJ 



1 uj 2 N 2 



(C.106) 



(C.107) 



4 4 

The spectrum for J 7 * will be slightly different, since complex conjugation is equivalent to 
taking n —> —n and q — > —q (or N — > —N, & b — >• — In this case, 



OJ 

2 



n+ (N - $>)e_ 



1 



UJ 



n- (N - $)e 6 - + 



a; 



(C.108) 



As is done previously, we have taken y, y*, z, and z* to be positive-definite, so that the 
corresponding fields do not diverge near the poles at x = ±1. 

Before computing the mode sum, it is helpful to introduce a notion of helicity similar 
to that used for the spin-1/2 field. If we define a quantity v such that £ = +l(£ = — 1) 
corresponds to the case of J- (J 7 *), then we can efficiently write 



3 + yt + zt + - 



1 uj 2 N' 
4 



(C.109) 



where 



UJ 



n + (N - $)e_ fe - b$ b + 



UJ 



UJ 



n-(N - <S>)e b - b$ b 



UJ 



(C.110) 



The gauge field is then understood to be comprised of two spin states with £ = 0, and two 
spin states with each of £ = ±1, respectively. 
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Mode Sum 

We wish to compute the small-i limit of the sum 



,% (u>, N, $ 6 , t) := S^(u, iV ' ^ *) = E E e_ 

5 £ J. 



(C.lll) 



Here, the additonal sum is over spin states, £ G {0,0,1,-1}. Given the spectrum derived 
previously, we see that two of the spectra will be identical to that of a scalar, and two will be 
slightly different. Similar what was done for the spin-1/2 field, we can make the identifications 



CO 



(or := § b - £<£[f where <3?jf = 1/uj) so that 



(C.112) 



(C.113) 



where A/" g f£ := oj(N g ^ — $ g f^) = M + 2^. In this case, S* g f can be written largely in terms of 
the scalar result: 

S gi (oj, N, t) = 2x S s (oj, N, $ 6 , t) + £ e *(A^ 2 -^)/4 5g(W) ^ ^ _ ^ t) 

£e{i-i} 

= 4 x S s ( w , TV, $ 6 , t) + ASgf (w, iV, t) . (C.114) 

In the last line above, we introduce the quantity A5 g f = ^^(S 1 ^ — S s ), which is a convenient 
one since the small-t limit of S^ f is very close to that of S s . (In what follows, we will dispense 
with the level of detail as was given for the spin-1/2 field and skip to the complete result, 
summed over spins.) Expanding out the argument of the sum above, we find a result that 
behaves similarly to the one found for the spin-1/2 field, in that the result will depend on 
whether the fluxes are either greater or smaller in magnitude when compared to <I>Q f = 1/u. 
When both \<$> b \ < <£q , we find that the differences, Asf , are given by As^ = 0, 



to 



-2oj + u 2 ^2\$ b \ 



Asf(u,N,<S> b ) = - 

OJ 

Asf(oj,N,<S> b ) = - 



(C.115) 
,(C116) 



ojjV 2 



+ oj z 



6 



OJ 

To 



. (C.117) 
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These expressions serve to cancel all the terms in s| which are odd under |<J>&| — > — |$&|, for 
those spin states with £ = ±1. 

Since the condition |$^| > <£g f = 1/co does not include very much parameter space for 
small u (since we are already implicitly assuming |3>;,| < 1), we will not consider this case in 
detail here. 

Spectrum and mode sum for massive gauge fields 

Let us now turn to the computation of the spectrum A™f f and the corresponding small-t limit 
of the mode sum 

Smgf(t) :=J>-' A " f (C.118) 

for a massive 6D spin-1 field A M on the rugby ball. 

The simplest field content for a massive gauge field includes, in addition to the metric and 
the gauge field, the scalar field <3? whose vev breaks the gauge group down to some subgroup 
H. The relevant part of the action is 

j d 6 x ^F MN F MN - (D M $? - U($)\ , (C.119) 

where the potential U($) is assumed to have a minimum for $ ^ 0. 

We would like to compute the linearized theory around a certain background solution 
with $ ^ (to give a bulk mass to some of the gauge fields) and the metric and gauge field 
set to an unwarped configuration with 4D Poincare invariance. So we substitute 

A M -> A M + V M and $ -> $ + rj , 

where V M and rj are small perturbations and now g M N, A M and $ represent the given back- 
ground solution. In order to interpret $ ^ as a 6D spontaneous symmetry breaking, we 
require to be constant and to lie at the minimum of U. Then, to solve the background 
scalar equation, D M D M $ = 0, we also demand that $ does not break the U(l) where the 
background gauge field lies: otherwise it would not be possible to set $ to a constant value, 



at least in the sphere compactification of interest in this paper [4C]. Such requirement is 
equivalent to demanding that the background gauge field lives in the Lie algebra of H. 

The linearized theory for the perturbations is rather complicated with mixing between 
V M , rj and the metric fluctuations. However, choosing the light-cone gauge defined in the 
previous subsection, the bilinear action for V M (such that Tr(A M V M ) = 0) and r\ has a 
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relatively simple form [37], that is 
d 6 x 



-g 



-D M V t D M V l - -D M V m D M V m 



1 rymn 

2 R 



VmVn 



&{T a ,T b }$ - gF mn V m x V n 

. d 2 u 



(C.120) 



i , d 2 u 



<9<I> 



:»7 



where the indices are raised and lowered by the background metric and e is a collective name 
for the gauge couplings of the generators broken by 

The simplest model of this sort having a non-trivial background flux is a £7(1) i x £7(1)2 
gauge theory with $ charged under the J7 (1)2 on ly, the background flux embedded in the 
£7(1) 1 and the potential having a simple mexican hat form, 



V|$| 2 + A 4 |$| 4 



(C.121) 



with /i 2 > (and so tachyonic) and A4 a positive constant. Let us first consider this simple 
option. The bilinear action for and 77 := (771 + ir\%)j\/2 is 



d 6 x- 



-D M V 2 D M V i2 



D M V 2 D M V m2 



e 2 $ 2 ,y2yi2 + V 2ym2 



R mny2y2 



^D MVl D M m - ^D MV2 D M m + m t 2 ach r? 2 - e 2 $ 2 r/ 2 



(C.122) 



where we have assumed <I> real and positive without loss of generality. The fields Vh and 
r/2 form together a massive 6D vector field with bulk mass m = 2|e|$, while 771 is a genuine 
scalar with bulk mass squared 2/i 2 . 

In the example above the massive gauge field is not charged under the background U(l), 
that is N = 0. A slightly more complicated non-abelian model can be used to illustrate the 
charged case. Let us consider for example an electroweak-like SU(2) x £7(1);, gauge theory, 
where <E> is in the 2i/ 2 representation. Like in the Standard Model we trigger SU(2) x £7 (1)& — > 
U(l) q through a non vanishing VEV, <£ T = (0, u/\/2), where U(l) q is the analogue of the 
electromagnetic £7(1), and we assume v real and positive without loss of generality. Also we 
embed the background gauge field in the U(l) q so that we can solve the scalar bulk equations 
for v constant and equal to the point of minimum of £7. In this case 



-±=(v + a + ir] Q ) 



a" 



V M = W£ I — + B M b, 



(C.123) 
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where a a are the Pauli matrices, the rji and a are real scalars and Wm an d Bm are the 
gauge field perturbations corresponding to the SU(2) and U(l)b group factors. By using the 
bilinear action in the light-cone gauge, eq. ( C.120Q , we find that = (W 1 ± iW 2 )/yJ~2 and 



77 = (771 ± iri2)/y/2 have the same bulk mass, m = \v\g2\, where 52 is the gauge coupling of 
SU(2), and represent altogether a massive gauge field with TV = ±1. Apart from this bulk 
mass term the light-cone-gauge bilinear action for W and r]^ coincide with that of a massless 
bulk gauge field and a massless scalar respectively. 

The bottom line of the examples above is that a massive gauge field leads to the 4D 
spectrum of a massless gauge field, given in eq. ( |4.5S| ), plus that of a scalar, which we 
provided before. It follows that the Sj coefficient of a massive gauge field are 

= 5s\ + Asf , s™ gf = 5s s 2 + Asf , (C.124) 

where the sf are the corresponding quantities for a 6D scalar, those given in eqs. d4.4j) — (|4.7|) , 
and we used \N\ < 1, to ensure the stability (see the appendix on the massless gauge fields). 
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